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We improve our description of tttt scattering data by imposing additional requirements to our 
previous fits, in the form of once-subtracted Roy-like equations, while extending our analysis up to 
1100 MeV. We provide simple and ready to use parametrizations of the ampUtude. In addition, 
we present a detailed description and derivation of these once-subtracted dispersion relations that, 
in the 450 to 1100 MeV region, provide an additional constraint which is much stronger than our 
previous requirements of Forward Dispersion Relations and standard Roy equations. The ensuing 
constrained amplitudes describe the existing data with rather small uncertainties in the whole region 
from threshold up to 1100 MeV, while satisfying very stringent dispersive constraints. For the 
SO wave, this requires an improved matching of the low and high energy parametrizations. Also 
for this wave we have considered the latest low energy Kt4, decay results, including their isospin 
violation correction, and we have removed some controversial data points. These changes on the data 
translate into better determinations of threshold and subthreshold parameters which remove almost 
all disagreement with previous Chiral Perturbation Theory and Roy equation calculations below 
800 MeV. Finally, our results favor the dip structure of the SO inelasticity around the controversial 
1000 MeV region. 
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I. INTRODUCTION 

In a series of papers [ll-ll| that we will denote by PY05, 
KPY06 and KPY08, respectively, we have provided sev- 
eral sets of precise phenomenological fits to tttt scattering 
data. The interest in a precise and model independent 
description of the data available in this process is twofold: 
On the one hand, it could be used at low energies to ex- 
tract information about the parameters of Chiral Pertur- 
bation Theory (ChPT) [4j] , quark masses and the size of 
the chiral condensate, pionic atom decays or CP violation 
in the kaonic system. On the other hand, in the interme- 
diate energy region, it could provide model independent 
information to identify the properties of hadronic reso- 
nances, particularly the scalar ones which are related to 
the spontaneous chiral symmetry breaking of QCD and 
the possible existence of glueball states. 

Pion-pion scattering is very special due to the strong 
constraints from isospin, crossing and chiral symmetries, 
but mostly from analyticity. The latter allows for a very 
rigorous dispersive integral formalism that relates the 
amplitude at any energy with an integral over the whole 
energy range, increasing the precision and providing in- 
formation on the amplitude even at energies where data 
are poor. Our aim is to provide reliable and model in- 
dependent TTTT scattering amplitudes that describe data 
and are consistent, within uncertainties, with dispersion 
relations. Note that, since we would like to test ChPT, 
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we are not using it in our analysis, and that, in order to 
calculate dispersive integrals up to infinity, we have been 
using Regge parametrizations obtained from a fit to data 
on nucleon-nucleon, meson-nucleon and pion-pion total 
cross sections In this work we will further improve 
our data analysis by imposing in the fits an additional 
set of once-subtracted dispersion relations, that we will 
also derive and describe in detail, showing that they are 
much more precise in the intermediate energy region than 
those we have used up to now. 

In general, on each paper of this series (or also in [6|) , 
we have first obtained a set of phenomenological "Uncon- 
strained" Fits to Data (UFD), which was fairly consistent 
with the dispersive requirements. Next, starting from 
that UFD set, we obtained "Constrained" Fits to Data 
(CFD) by imposing simultaneous fulfillment of dispersion 
relations. These constrained fits not only describe data, 
but are remarkably consistent with the strong analyticity 
requirements. Furthermore, the output of the dispersive 
integrals is model independent and very precise. 

The constraints we imposed in the first two papers of 
this series were just a complete set of Forward Dispersion 
Relations (FDR), plus some crossing sum rules. In the 
third paper, apart from including the most recent and 
reliable data up to that date on Ki4, decays 0, Hi, we 
also imposed Roy Equations [9*], because they constrain 
the t ^ Q behavior of the amplitude, while ensuring s — t 
crossing symmetry. These equations, which had already 
been used in the 70's to analyze some of the existing data 
[To| , as derived by S. M. Roy, have two subtractions and 
provide a strong constraint in the low energy part of the 
partial waves. For this reason there has been recently a 
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considerable effort to analyze them in relation with ChPT 
They have also been recently used to eliminate [12] 
the longstanding ambiguity about "up" or "down" type 
solutions of the SO wave data analyses. Since Roy equa- 
tions are written in terms of partial waves, they lead, if 
supplemented with further theoretical input from ChPT 
[13| , to precise predictions for resonance poles like the 
much debated /o(600). Despite being listed with huge 
uncertainties in the Particle Data Book [3] , several anal- 
yses using analytic methods or dispersive techniques with 
chiral constraints [1, as well as those using Roy Eqs. 

, are in fair agreement about its pole position, around 
450 — z250 MeV. However, its nature remains controver- 
sial, since it might not be an ordinary meson [l6l |. A 
precise analysis of tttt scattering data may help clarifying 
the situation by studying the /o(600) parameters (like the 
coupling lyWi- and the connection of the pole to QCD pa- 
rameters fla] , although one have to bear in mind Jl9j the 
difficulties to interpret the coupling in terms of simple 
intuitive models. Nevertheless, let us remark that here 
we only aim at a precise description of data, which could 
later be used for those purposes among many others, but 
the interpretation of this resonance and the extension to 
the complex plane are beyond the scope of this work. 

Back to Roy equations, when used only with data, as 
it is our case, the 82 wave scattering length, which is 
very poorly known experimentally, dominates completely 
the Roy equations uncertainties, that become very large 
above roughly 450 MeV, for the SO and S2 waves. For 
that reason Roy equations do not provide a significant 
additional constraint for the amplitudes beyond that en- 
ergy, once they are already constrained with FDR. In this 
work we will overcome that caveat with additional once- 
subtracted Roy-like equations that have a much weaker 
dependence on scattering lengths. The fact that these ad- 
ditional equations have a much smaller uncertainty above 
roughly 450 MeV will force us to refine the matching of 
our SO wave parametrizations. 

Let us remark, though, that our parametrizations are 
consistent with those in KPY08 within one standard de- 
viation, with the only exception of the SO wave. However, 
the new central values satisfy Roy equations and the new 
once subtracted dispersion relations better. Moreover, 
we will now be able to extend the Roy equations analy- 
sis, both with one and two subtractions, up to 1115 MeV, 
instead of just the KK threshold. 

Once again we remark that the functional form of the 
amplitude parametrizations becomes irrelevant once the 
imaginary part of the amplitude is used in the disper- 
sive integrals, whose results are model independent. In 
the understanding that running the dispersive represen- 
tation could be tedious for the reader, we provide results 
in terms of our simple and ready to use CFD parametriza- 
tions, which are very good approximations to the disper- 
sive result. 

The plan of this work goes as follows: In Section [H] 
we very briefly comment on the simple unconstrained 
data fits (detailed in Appendix A) of all partial waves 



obtained in previous works. Only the SO wave is given in 
more detail in Section [IIII to introduce the new improve- 
ments. These are of two kinds: On the one hand, the 
data has changed, since we are taking into account the 
final and more precise NA48/2 data [20], including the 
threshold enhanced isospin violation correction to all K^i 
data, and getting rid of the controversial K 2tt datum. 
On the other hand, we have improved our parametriza- 
tion, by imposing a continuous derivative matching be- 
tween the low and intermediate energy regions and al- 
lowing for more flexibility in the parametrization around 
the /o(980) region. 

In Section ITVl after introducing FDRs and Roy equa- 
tions very briefly, we present the once-subtracted dis- 
persion relations and compare their structure with the 
standard Roy equations. Next, in section |V] we impose 
these new relations together with the constraints already 
used in previous works (FDRs, sum rules, standard Roy 
equations...) to obtain the final representation for the 
amplitudes, i.e., the CFD set of amplitudes. In Sect I VII 
we study the threshold parameters and Adler zero de- 
terminations stemming from this constrained fit through 
the use of additional sum rules and dispersive integrals. 
Then, in the discussion section, we compare these CFD 
with our previous results and other works in the litera- 
ture, and we comment on the effect of considering differ- 
ent choices of data or parametrizations as a starting point 
to obtain our final result. In particular, we show how our 
results favor a "dip" structure in the SO wave inelastic- 
ity right above 1000 MeV, which has been the subject of 
a longstanding controversy [2l|. Finally, we present our 
conclusions. In the Appendices we provide a list of all 
parametrizations and parameters of the UFD and CFD, 
as well as the detailed derivation of the once-subtracted 
relations together with all relevant integral kernels. In 
appendix |D] we provide a table with the phase shifts in 
the elastic region, as obtained from the dispersive repre- 
sentation. 



II. THE UNCONSTRAINED FITS TO DATA 
A. Our previous works 

To explain the motivation for further improvements 
in our previous amplitudes, we briefly describe next the 
results of the previous articles. In particular, 

- In PY05 [1] we obtained simple and easy to use 
phenomenological parametrizations of tttt scattering data 
whose consistency was checked by means of FDR and 
several crossing sum rules. The P, S2, DO, D2, F, GO 
and G2 partial waves were described by simple flts to tttt 
scattering data up to 1.42 GeV. In the elastic regime, 
the P wave was obtained from a fit to the pion form 
factor. For the SO wave, given the fact that there are 
several conflicting sets of data, we fltted first each set 
separately and then performed another global fit only in 
the energy regions where different data sets are consis- 
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tent. Surprisingly, some of the most commonly used data 
sets failed to pass these consistency tests, although the 
global fit was in fairly good agreement with FDR. Hence, 
it could be used as a starting point for a constrained fit 
to data. This CFD was obtained by imposing FDR and 
crossing sum rules to be satisfied within errors, in the 
elastic regime and up to 925 MeV. As a result, a precise 
description of the data up to 925 MeV was obtained by 
means of a constrained fit, satisfying the FDR and sum 
rule requirements remarkably well. 

- In KPY06 [21 we refined our parametrizations above 
KK threshold, including more tttt data but, most impor- 
tantly, TTTT — T' KK data in a coupled channel fit. These 
reduced uncertainties forced us to slightly refine the UFD 
parametrizations of our DO, D2 and P waves between 1 
and 1.42 GeV as well as the Regge parameters. This led 
to a remarkable improvement in the consistency of the 
7r°7r" FDR. 

- In KPY08 [1] we also considered Roy equations ^] 
for our amplitudes below KK threshold. The UFD fits, 
where we had previously incorporated ^ the most reli- 
able low energy data from Kn decays to that date Q, 
satisfied Roy equations fairly well and the agreement was 
remarkably good once they were imposed into a new set 
of CFD. 

Since, in this work, we are going to consider a set of dis- 
persion relations in addition to the dispersive constraints 
we have just described, our starting point will be the 
UFD set already obtained in KPY08, that we describe 
only very briefly in the next subsections, but explain in 
detail in Appendix A. The only exception will be the SO 
wave, that we describe in Sect. IIIII The reasons are the 
appearance of new data (20| . the existence of modifica- 
tions on the analysis of the old experimental results, and, 
in addition, that we have found that the new constraints 
are strong enough to require a better matching, with a 
continuous derivative, between the low and intermediate 
energy parametrizations. 



B. Notation 

For TTTT — > TTTT Scattering amplitudes of definite isospin 
/ in the s-channel, we write a partial wave decomposition 
as follows: 



TT ^ — ' 



(1) 



where S'^p (s) and (s) are the phase shift and inelastic- 
ity of the /, (. partial wave, i is the angular momentum, 
and k is the center of mass momentum. In the elastic 
case, Tj — \ and 



/f'(.) = sin4^)(.)e<'(^). 



Note that 7 = 0,1,2 and that whenever I is even (odd) 
then i is even (odd), and thus we will omit the isospin in- 
dex for odd waves. We may refer to partial waves either 
by their /, £ quantum numbers or by the usual spectro- 
scopic notation SO, S2, P, DO, D2, F, GO, G2, etc... 

In addition, we recall the expressions for the so called 
threshold parameters, which are the coefficients of the 
amplitude expansion in powers of center of mass (CM) 
momenta around threshold: 



,1/2 



(I) 



dl) 1.2 



2M, 



(3) 



Note that and are the usual scattering lengths 
and slope parameters. Customarily these are given in 
A/tt units. 



C. Parametrizations for S2, P, D, F and G waves 

The S2, P, DO, D2, F and G waves are described by 
very simple expressions. For the S2, P and DO waves 
we use separate parametrizations for the "low energy 
region", i.e. energies s^/^ < s]^"^ ~ IGeV, and the 
"intermediate energy region" , which extends from the 



1 /2 

matching energy Sj^^ up to 1.42 GeV. For each wave, 

s^j"^ is typically the energy where inelastic processes can- 
not be neglected. Note that, above 1.42 GeV we will 
assume that tttt amplitudes are given by Regge formulas, 
which correspond to fits to experimental data (see [Sj] and 
KPY06 for details). 

In the "low energy region" , where the elastic approxi- 
mation is valid, we use a model independent parametriza- 
tion for each partial wave i^^', that ensures elastic uni- 
tarity: 



t 



(I) 



1 



2k 



cot 6^ (s) — i 



To ensure maximal analyticity in the complex plane 
coti5^^'*(s) is then expanded in powers of the conformal 
variable 



w{s) 



(2) 



where Si is a convenient scale for each wave, to be pre- 
cised later, always larger than the s range where confor- 
mal mapping is used. The use of a conformal variable al- 
lows for a very rapid convergence — at most two or three 
terms are needed in the expansion — so that each wave 
is represented by only 3 to 5 parameters, corresponding 
to the coefficients of the expansion and the position of 
the zeros and poles when we have found convenient to 
factorize them explicitly @. We remark again that the 
use of a conformal expansion does not imply any model 
dependence. 

In the intermediate energy inelastic region, we have 
used purely polynomial expansions both for the phase 
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shifts and inelasticities in terms of the typical energy or 
momenta involved in the process. 

All these simple parametrizations have been fitted to a 
large number of experimental data on tttt phase shifts or, 
in the case of the P wave, to the vector form factor data, 
which gives much more precise results. In Appendix A, 
we provide the detailed parametrizations for each partial 
wave, together with the resulting parameters and their 
uncertainties, from now on denoted by p^^'' and Spi, re- 
spectively. 

Let us remark that, on a first step, each partial wave 
has been fitted independently of each other, without im- 
posing any constraint from dispersion relations, and that 
is why we refer to such initial fits as "Unconstrained" 
Fits to Data or UFD. In KPY08 we showed that these 
UFD provided a good description of data, and a fairly 
reasonable consistency in terms of dispersion relations. 
Of course, the consistency is much better, remarkable 
indeed, once we impose the dispersion relations as con- 
straints to the fit, but then all waves become correlated. 
The uncorrelated fits, apart from providing the starting 
point of our calculation, and although they are less re- 
liable than our final constrained results, could be of rel- 
evance if new and more precise data becomes available 
for a given partial wave, since then only that particular 
partial wave should be modified, without affecting the 
others. 



III. SO WAVE PARAMETRIZATION 

This is the only wave that changes in the new sets of 
unconstrained data fits. This is due to three reasons that 
we will explain in separate subsections. 



A. On isospin violation in K14 decays 

There has been a recent calculation [2^ showing that, 
due to threshold enhancements, isospin corrections in 
Ki4 decays 0, H, [131 could be larger than naively ex- 
pected. A leading order ChPT calculation has been pro- 
vided to correct the phase shift determination in the 
isospin limit, that should be valid within the whole range 
of Ki4 decays. Note that the uncertainties in the previ- 
ous UFD set in ^ were obtained taking into account 
systematic errors on the data, including possible isospin 
corrections, but only of natural size. Since the most re- 
cent data from Ku decays play a relevant role in the SO 
wave of our UFD set, and the suggested isospin breaking 
effect is unnaturally large, we will modify the SO wave 
by correcting the K14 data as suggested in [l^l, so that 
it can be used in our isospin limit formalism. Note that 
this isospin correction was already made available in Q 
and again in the final NA48/2 results [13]. 



B. The K ^2n data 

Let us emphasize again that this is a data analysis, 
and, as such, it depends on whether we include or not 
certain experimental results that are somewhat contro- 
versial. This is for instance the case of the phase shift 
difference obtained from K ^ 2tt decay [l^ that we used 
in KPY08: 

(57.27 ± 0.82exp. ± 3rad. ± lchPTappr.)° • (4) 

The extraction of the tttt scattering phase from this de- 
cay is affected by large uncertainties that have to be es- 
timated from ChPT. A similar value is obtained if using 
the Particle Data Group data and the prescription for 
radiative corrections in In Q we took the simple 

linear sum of the errors quoted in (23j . which is larger 
than the usual quadrature addition. However, the use of 
the datum above has been questioned in [25|, also sug- 
gesting that it could be partly responsible for the differ- 
ences between our approaches in the intermediate energy 
region. It is true that this data point always lies some- 
what above our parametrizations of KPY08, 51.7 ± 1.2° 
for the UFD and 50.4 ± 1.1° for the CFD, and even more 
so from those in [ll|, 47.7 ± 1.5°. While preparing this 
work, a re-analysis has appeared [26| taking into account 
more precise experimental data and other improvements 
including an update of the low-energy constants, yield- 
ing: 

SqHm],) - (5(')(M|) - (52.5 ± 0.8exp. ± 2.8theor.)°. (5) 

This is still compatible with the value in Eq. Q, but 
seems in much better agreement with tttt scattering de- 
terminations. However, this new extraction uses as an 
input the SO phase shift value from a tttt scattering analy- 
sis usiiig Roy equations and ChPT, obtained by the Bern 
group [llj . Thus it would be somewhat circular to use it 
as input in our approach. Furthermore, we have studied 
the alternative scenarios with and without the K ^ 2t: 
value in our fits, finding that the scenario without it is 
slightly preferred by dispersion relations. For these rea- 
sons, we will present results for fits removing the — >■ 27r 
controversial datum. As a consequence, our new uncon- 
strained fits have somewhat smaller errors than those in 
KPY08, which makes dispersion relations harder to be 
satisfied. 



C. Improved parametrization and matching 
condition between low and intermediate energies 

In previous works, only continuity, but not a contin- 
uous derivative, was imposed for the SO phase shift at 

1/2 

the matching point, then chosen at s/j = 932 MeV . 
It has been suggested [s^ that such a crude matching 
could explain the roughly 2a level discrepancies in the 
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SO wave between the KPY08 analysis and that of the 
Bern group 
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in the 450 - 800 MeV region. We have 
checked that the improved matching by itself only af- 
fects the SO wave sizably in the /o(980) region, although 
the effect is rather small below. However, this improved 
matching adds together with the two effects in kaon de- 
cays discussed above, to become a relatively larger effect 
that certainly improves the agreement with the predicted 
SO wave in 

In this work we want to keep the same low energy con- 
formal parametrization of KPY08 or j^]. However, to im- 
prove the flexibility of the parametrization we will keep 
one more term in this expansion. Actually, it has been 
pointed out that the difference between the parametriza- 
tion in KPY08 and that of ll'] could be due to the fact 
that our conformal parametrization at low energies was 
not sufficiently flexible [i^. The additional parameter 
does not improve significantly the fulfillment of disper- 
sion relations nor the data fit, but the output of the dis- 
persion relations with one parameter less would violate 
very slightly the elastic unitarity condition around 500 
MeV. For that reason we keep this additional term, and 



cot 6^°\s) 



2k s - ^Zq 



(6) 



— ^ + Bo + B^w{s) + B2w{sf + B^w{sf 



^^)- 777^ , .o = 4A4. (7) 

where the new values for the UFD parameters are: 

Bo = 7.26 ± 0.23, Bi = -25.3 ± 0.5, 
B2 = -33.1 ± 1.2, B3 = -26.6 ± 2.3, zq = (8) 

which are obtained with the same procedure as in (J| but 
now including the additional -B3, the isospin corrections 
and getting rid of the K ^ 27r data, as already com- 
mented in subsections IIII Al and IIII BI above. Namely, 
in this fit we have considered the data on K14 decays 
[Ti], including the final Ki4 decay data from NA48/2 [23| 
(which supersedes Q), and a selection of all the existing 
and often conflicting tttt scattering data [27", '2^ . This 
selection corresponds to an average of the different ex- 
perimental solutions that passed a consistency test with 
Forward Dispersion Relations and other sum rules in the 
initial work PY05. To this average we assigned a large 
uncertainty to cover the difference between the initial 
data sets. For the sake of brevity we simply refer to that 
work, or the Appendix of Ref. @, for a complete and 
detailed description of the data selection. Uncertainties 
in Eq.® come from data only. In order to use the UFD 
by itself, a systematic uncertainty due to parametriza- 
tion dependence [29"] should be taken into account. But 
as we have seen, possible parametrizations are strongly 



restricted by imposing dispersion relations and unitarity 
in their output, thus reducing dramatically this source 
of systematic uncertanty. Hence, we will only quote the 
data uncertainty for the CFD. Of course, since dispersion 
relations are imposed within uncertainties, the residual 
parametrization dependence is reflected in the error bars 
from the result of the dispersive representation, which we 
give in Table XII of Appendix D. 

Despite this amplitude being used only in the physi- 
cal region, we have explicitly factorized a zero at sa = 
= Mll2 ~ (98.7 MeV )2 for these unconstrained 
fits. This corresponds to the position of the so called 
Adler zero, required by chiral symmetry [30], at lead- 
ing order in ChPT. Note, however, that this zero lies 
very close to the border of the convergence region of the 
conformal expansion (see Fig. 16 in KPY08), which is 
therefore not very well described by the expression above. 
Hence, zq should not really be interpreted as the exact 
position of the Adler zero, but just as another parame- 
ter of our parametrization. Of course, the physical low 
energy region, which is the only one relevant for the dis- 
persive representation, lies well inside the convergence 
region of the conformal expansion, and is very well de- 
scribed by Eq. ([5]). Actually, we will show in Sect. IVII 
below that, when this parametrization is used inside the 
dispersive representation, one finds an Adler zero in the 
correct position. 

Let us now turn to the intermediate energy region. In 
previous works, a two-channel K-matrix formalism, fol- 
lowing the experimental reference in j28j , was used to de- 
scribe the region around KK threshold. This is a rather 
popular formalism to describe multichannel scattering of 
two-body states, but has several disadvantages for our 
purposes. One, of course, is the use of only two channels 
TTTT and KK, neglecting possible inelasticity contribu- 
tions from 47r or other channels. These are rather small, 
but since we aim at a precision determination, we should 
allow for more flexibility on the inelasticity, whereas the 
two channel K-matrix yields a strong relation between 
phase and inelasticity. The second caveat is the huge 
correlations between K-matrix parameters, which makes 
it very hard to improve by means of constrained fits, as 
we will do later on. Finally, a very strong disadvantage is 
that the phase dependence on the K-matrix parameters 
is so complicated that it is not possible to make an ana- 
lytic matching with the low energy parametrization, and 
a numerical matching is much more ineffective and harder 
to implement. Let us note that some of these caveats 
were already removed when using some very naive poly- 
nomial parametrizations considered in the Appendix B of 
KPY06. We will use those same parametrizations here 
but with additional terms in the expansion to compen- 
sate the loss of flexibility due to the improved matching 
conditions. In particular, between the matching point 
and 1.42 GeV, we will use: 
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do 



M 



B 



K 



Ml 



Ml 



+ o 



Ml 



K 



where k2 = - Ml, = ^s/i-M^ and do 

is the phase shift at the two kaon threshold. Note, 
however, that we have lowered the matching point to 

1/2 

s^.j — 850 MeV, since we have found empirically that 
this helps improving dispersion relation ftilfillnicnt, as the 
slope is somewhat smaller there. As a final remark, we 
have added a term proportional to the rj momentum, to 
reflect the opening of the rjr] channel, which has shown to 
have some relevance in the description of the data [3l| . 
In this respect we want to clarify a common source of 
confusion about Roy (or GKPY) equations: These rela- 
tions include all possible coupled channels contributions, 



(0.85GeV)2 <s< 4Ml, 

(9) 

4M2)-^, ami < s < (1.42GeV)2 

I 

or at least are consistent with them, as long as they are 
in agreement with the experimental inelasticicity. This 
simple term is purely phenomenological, and given the 
size of the experimental errors this additional term is 
more than enough to just describe the cusp due to the 
presence of this channel. However, it yields very slightly, 
but favorable, difference in the fulfillment of dispersion 
relations. 

By defining Sm = S{sm) and (5^/ = dS{sM)/ds, which 
are obtained from Eq. ([S]), and kM — 1^2 (sm)!, it is rather 
straightforward to impose continuity and a continuous 
derivative for the phase shift at sm, to find: 



\k2_ 
kM 



do{l-T^] +Sm^{2 



M 



K 



CM 



rr "-2 



^) + |fe|(fcM 



^3 



— — + D6(s - 4M ) 



K 



^(fcMjM^) ^ (0.85 GeV)2 < s < AMI, 



Ml 



AMI < « < (l-42GcV)2. 

m 



As previously commented, with the exception of the 
K ^ 2iT datum, the inclusion of isospin corrections 
to data explained above, and our use of the final 
NA48/2 results (l^, our treatment and selection of data 
for the phase is exactly the same one followed in the 
previous works and so we will not repeat them 
here. In Table |V] of the Appendix, we provide the val- 
ues for the do,c,B,C and D parameters resulting from 
the unconstrained fit to those data. In Fig. [T] we show 
the resulting phase from the unconstrained data fit to 
the SO wave phase shift up to 1420 MeV, and in Fig. [2] 
we show the low energy re gion in detail, including the 
isospin violation correction ^22] that we have subtracted 
from all the Kn data. Note that this correction amounts 
to slightly less than 1 degree in the region from threshold 
to 400 MeV, which is not much at high energies, but very 
relevant close to threshold. 

In Fig. [3] we show a comparison of the phase shift re- 
sulting from the new UFD with the improved matching 
versus the one obtained in KPY08. The changes at low 
energy are due to the update on the Ki4^ data and their 
isospin corrections, together with the fact that we now 
discard the K ^ 271 datum. The bump in the 500 to 



800 MeV region observed in KPY08 has almost disap- 
peared. Thus, the improvement on the data and its cor- 
rections reduces almost completely the disagreement of 
our UFD description with the phases in [ll| — the line la- 
beled CGL in the plot — although our central values are 
still larger in the 550-800 MeV region. Furthermore, as 
we will see later, for the constrained fits we are in an 
even better agreement with (TJ]. The changes above 
the matching point are sizable for the phase, mostly 
around the sharp phase increase usually associated with 
the /o(980) resonances, as can be seen in Fig. [3] where 
the central value for the new phase is compared with 
that in KPY08. Note the much smoother behavior in 
the matching region for the new UFD parametrization 
and the more dramatic KK threshold effect. 

Concerning the SO wave inelasticity, we approximate it 
to 1 up to the two-kaon threshold, and use the following 
parametrization above that energy: 



exp 
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FIG. 1: The new SO wave unconstrained fit (UFD), where the dark band covers the uncertainties, versus the existing phase 
shift data from [27l. [2^. Note that the if — > 27r point has been excluded from the fit as explained in the text. 
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FIG. 2: The new SO wave unconstrained fit (UFD), where the 
dark band covers the uncertainties, versus the "old" phase 
shift data from Ku decays 0] together with the final NA48/2 
results, which supersedes the data from the same experi- 
ment [91 that we used in KPY08. We are also showing the 
isospin violation correction [23|, which has been included in 
the data shown here. Finally, we show the results of the CFD 
parametrization to be explained in Sect. [V] which is almost 
indistinguishable from the UFD curve. 



for 4M|. < s < (1.42 GeV)^. By neglecting the term pro- 
portional to the 77 momentum, which is numerically very 
small as seen in the Appendix A.l, and by re-expanding 
the above equation in powers of k2/s^^'^ up to third order, 



we recover the polynomial expression in KPY06, but the 
definition above ensures the < 77q°'' < 1 physical condi- 
tion, whereas the simple polynomial in KPY06 did not. 

For the inelasticity data, we follow again the same se- 
lection as in previous works of this series, but we do not 
include now the data from Kaminski et al. [l^] in the 
calculation; we only consider the 1973 data of Hyams 
et al. [13] and Protopopescu et al. 27 1. The reason is 
that the main source of uncertainty is systematic, and 
if we include the large number of points of Kaminski et 
al. with their huge statistical errors, the outcome of the 
fit has much smaller errors than the original systematic 
uncertainties. By keeping only the other two sets, which 
are incompatible, we obtain a fit with a large /d.o.f., 
and by rescaling the uncertainties in the inelasticity pa- 
rameters we mimic the dominant systematic uncertain- 
ties much better. Of course, our results are still in very 
good agreement with Kaminski et al. Was the systematic 
uncertainty not dominant, this would not be necessary. 
In Table |V] of the Appendix, we provide the values for 
the ii parameters, and in Fig. |4]we show the results of 
the unconstrained fit to the SO wave inelasticity data up 
to 1420 MeV. 

Finally, let us remark that the inelasticity is the scat- 
tering parameter that suffers the biggest change with re- 
spect to the KPY08-KPY06 parametrization, as can be 
seen in Fig. [5l The new parametrization shows a big dip 
in the inelasticity between 1 and 1.1 GeV, whereas the 
KPY08 one does not. As already commented in PY05, 



this is a longstanding controversy (see, for instance 21 1 
and references therein) between different sets of data 
coming from pure tttt — > tttt scattering versus those com- 
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FIG. 3: Fit to the SO wave phase shift, with the improved continuous derivative matching (UFD, continuous hue) versus the 
simpler one used in KPY08. We also show the phase predicted in [ll[ (CGL). 
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FIG. 4: The new SO inelasticity fit (UFD set) to the tttt tttt 
scattering data of Hyams et al. (1973) and Protopopescu et 
al. As explained in the text, we do not fit the Kaminski et al. 
data [231 although our fit is compatible with them. The dark 
band covers our uncertainties. For all data sets see Fig. 1181 



ing from tttt — KK analysis. Actually, in PY05 (see 
Fig. 6 there) we considered both possibilities: we found 
that Forward Dispersion Relations favored the "non-dip 
solution" very slightly, but we kept the " dip-solution" in 
order to use the phase and inelasticity coming from the 
same experiment. In KPY06 we found a similar situation 
but since the K-matrix slightly preferred again the "non- 
dip solution", this time we decided to use it. However, 



in terms of fulfillment, the difference is minute for FDRs, 
and even more so for standard Roy equations, since, as we 
have already commented and we will see in detail below, 
the uncertainties in the subtraction constants become so 
large above 500 MeV that we cannot use them to discard 
any of the two scenarios. The existing set of dispersion 
relations did not allow us to make a really conclusive 
statement about the inelasticity in the 1 GeV region. 

One of the main results of this work is the deriva- 
tion and use of once subtracted Roy-like dispersion re- 
lations, the GKPY equations presented in Sect. lIVDl be- 
low, which are more precise in the 1 GeV region and 
clearly favor the solution with a dip, thus helping to set- 
tle this dip versus non-dip controversy. 



IV. DISPERSION RELATIONS AND SUM 
RULES 

From the theoretical side, tttt scattering is very spe- 
cial due to the strong constraints from isospin, crossing 
and chiral symmetries, but mostly from analyticity. The 
latter allows for a very rigorous dispersive integral for- 
malism that relates the tttt amplitude at any energy with 
an integral over the whole energy range, increasing pre- 
cision and providing information on the amplitude even 
at energies where data are poor, or in the complex plane. 

Let us emphasize once more that the dispersive ap- 
proach is model independent, since it makes the data 
parametrization irrelevant once it is included in the inte- 
gral. The previous works [1, @| of this series made use of 
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FIG. 5: Fit to the SO wave inelasticity (UFD) with the im- 
proved continuous derivative matching (continuous Hue) ver- 
sus the simpler one used in KPY08 (dashed line). The dark 
band covers the uncertainties of the former, whereas the dot- 
ted curves enclose the uncertainties of the latter. Note that 
the drop in the inelasticity right above 1 GeV has become 
much deeper. In contrast to Fig. [J] we only show the data 
coming from tttt — > KK and the tttt — >■ nn on which is based 
the KPY08 fit. For aU data sets see Fig. [H 



two complementary dispersive approaches, Forward Dis- 
persion Relations and Roy equations, that we briefly re- 
view next, before introducing the new set of once sub- 
tracted Roy-like equations. 



A. Forward Dispersion Relations (FDR) 

They are calculated at t = 0, so that the unknown 
large-t behavior of the amplitude is not needed. There 
are two symmetric and one antisymmetric isospin com- 
binations to cover the isospin basis. For further con- 
venience we will write them as a difference Ai{s) that 
should vanish if the dispersion relation is satisfied ex- 
actly. In particular, the two symmetric ones, for Tr^Tr"*" 
and 7r°7r°, have one subtraction and imply the vanishing 
of 



A,{s) = ReF,{s,0)- F,{AMlO)- 



s{s - 4Af2) 



RP. 



(2s'-4M2)ImF,(s',0)ds' 
4M? s'(s'-s)(s'-4M2)(s' + s-4M2) 



(12) 



where Fi stands for the Fo+(s,t) or Foo{s,t) amplitudes, 
and "P.P." stands for the principal part of the integral. 
They are very precise, since all the integrand contribu- 
tions are positive. The antisymmetric isospin combina- 
tion /t = 1 does not require subtractions and implies the 



vanishing of the following difference: 

2s - 4M2 



A(,,=i)(s)^ F^^-^Us,0)- 
P.P, 



(13) 



°° , , ImF(^*=i)(s',0) 
d s 



I^Mi (s'-s)(s' + s-4M2)- 
All FDRs are calculated up to = 1420 MeV. 

B. Roy Equations 

These are an infinite set of coupled equations [9] , equiv- 
alent to non-forward dispersion relations plus t — s cross- 
ing symmetry. They are well suited to study poles of 
resonances and scattering data, since they are written 
directly in terms of partial waves t^/^ of definite isospin 
/ and angular momentum £. Remarkably, S. M. Roy 
managed to rewrite the complicated left cut contribution 
as a series of integrals over the physical region. In the 
original work of Roy and all applications until now, the 
convergence of the integrals was ensured by making two 
subtractions. 

As we did with FDR, we will recast each one of the 
Roy Equations as the difference 



41) 



i'=o e'=o 



4M2 



ds'K',^,{s,s')lmtj,{s'), (14) 



that should vanish when the equation is exactly satisfied. 
Roy equations provide as output the real part of partial 
waves below 1115 MeV. Although, in principle, one could 
consider output for waves up to higher £, in this work we 
are interest in results for f = 0, 1 only. Hence, we have 
separated those waves explicitly below Smax- 

As it was done in KPY08, below sUax = 1420 MeV, we 
consider the imaginary parts from all our £ < A partial 
wave parametrizations as input. Above that energy, we 
take into account all waves together parametrized with 
Regge theory — see Appendix IA8I The K^l (s, s') are 
known kernels, and thus we will refer to the integral 
terms as "kernel terms" or KT{s). The "driving terms", 
DT^ (s), have the same structure as the kernel terms, but 
their input contains both the contribution from £ — 2,3 

1/2 

partial waves up to Smax = 1420 MeV, and the Regge 
parametrizations above. We have explicitly checked that 
the £ — A contribution below Smax is irrelevant, so that 
we will refer just to waves up to ^ = 3. Finally, the 
so-called subtraction terms are given by: 



4Afi 



12 A/2 



(15) 



(2a[! - 5al){SioSm + -^SnSi 



6 



It is very relevant to remark once more that these equa- 
tions have two subtractions, as can be seen by the pres- 
ence of the term proportional to (s — AM^){2aQ — Soq). 
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This strong energy dependence of ST{s) makes these 
twice subtracted Roy Equations very suitable for low en- 
ergy studies, and even more so when complemented with 
theoretical predictions of the scattering lengths coming 
from ChPT [ll|- 

Roy Equations are vahd up to -^s < SA/tt — 1120 MeV. 
However, we will see that the uncertainties in the scatter- 
ing lengths, when propagated to high energies, become 
too large above roughly 450 MeV, due to the term pro- 
portional to s. For this reason, in KPY08 it did not 
make sense to deal with the complications of a precise 
description around KK threshold and thus we imple- 
mented them up to 2Mk- One of the main novelties of 
the present work is that, since the once-subtracted Roy- 
like equations explained below will have much smaller un- 
certainties in the KK threshold region, we have now im- 
plemented these new equations, together with the stan- 
dard Roy equations, up to 1115 MeV. 



C. Two Sum Rules 

Apart from FDRs and Roy equations, two sum rules 
that relate high energy (Regge) parameters for i to 
low energy P and D waves, have been considered through- 
out previous works. 

The first sum rule (PY05) is nothing but the vanishing 
of the following difference 

/ - n ds Ii^^^^'"^Hs,4Af2) - ImF(^*=i)(s,0) 
8M^[s~2M^]lmF('^='\s,0) 

4m/^ s^is-iM^.y ' ^''^ 

where the contributions of the S waves cancel and only 
the P and D waves contribute (we also include F and 
G waves, but they are negligible). At high energy, the 
integrals are dominated by the rho reggeon exchange. 

The second sum rule we consider is given in Eqs. (B.6) 
and (B.7) of the second reference in |Tl|, which requires 
the vanishing of 



J= I ds 



■ 4ini (s, 0) - lOIm (s, 0) 

s2(s-4A/2)2 



_ ,^ ImF^(i)(5,0)-Imf(i)(s,0) 



(17) 



Here, F'^^'>{s,t) = dF'^^^s,t)/d cos 6. At high energy, 
the integral is dominated by isospin zero Regge trajecto- 
ries. 



D. GKPY Equations 

The main novelty of this work is that we present and 
use a new set of Roy-like dispersion relations for tttt 



scattering amplitudes. For brevity, we will call them 
GKPY equations, as we have already done when present- 
ing some partial and preliminary results in several confer- 
ences jsil, ^33^ . In brief, their derivation follows the same 
steps as for Roy equations, starting from fixed t disper- 
sion relations for a complete isospin basis, that S. M. Roy 
subtracted twice to ensure that the integrals converged 
when extended to infinity. However, by using the com- 
plete set of isospin amplitudes Fqq, Fo+ and it 
is easy to see that one subtraction is enough. Actually, 
the two first amplitudes are s — u symmetric and the 
contributions from the s and u channels, that would be 
divergent by themselves alone, cancel when considered 
simultaneously. The f (^*=i) amplitude is dominated by 
the rho Regge exchange and neither the left nor the right 
cut are divergent with one subtraction. We provide the 
detailed derivation in Appendix |B1 which leads to the 
vanishing of the following difference: 



^GKPY (/) _ ^(/) ^ST\-DT\ {s) 
-EEP-P-/ ds'K\l,{s\s)lTnt^:\s'). (18) 

The subtraction terms STg are linear combinations of 
scattering lengths Oq, and can be found in Appendix|B] A 
very relevant observation for this work is that, in contrast 
to the standard Roy Equations, the subtraction terms in 
GKPY do not depend on s. 

The integral and driving terms DT (s) in Eq. (jl8p are 
analogous to the kernel and driving terms in Roy equa- 
tions, but the integrals contain the k\1, kernels, instead 

of the if// . The explicit expressions for k\1, are lengthy 
and we provide them in Appendix [C] Note that, as the 
once subtracted GKPY equations have kernel terms that 
behave as ^ l/s2 at higher energies, instead of the ~ 1/s'^ 
behavior in Roy equations, the weight of the high energy 
region is larger. Nevertheless, the contribution to the 
driving terms coming from energies above 1.42 GeV is 
generically smaller than the contribution coming from 
the D and F waves below 1.42 GeV, which means that 
their infiuence is still under control. 



E. Roy versus GKPY Equations 

Fig. |6] presents a decomposition of Roy equations for 
the SO, P and S2 waves into four parts: the "in" part, that 
represents what our parametrizations give for Re t^P , the 
subtracting terms ST{s), the kernel terms KT{s) and the 
driving terms DT{s). Note that, for these equations to be 
satisfied exactly, the first contribution should equal the 
sum of the other three. The numerical calculations have 
been performed by taking the UFD amplitudes described 
in the previous Sections as input. For illustration, we 
have drawn as a gray area the region that violates the 
unitarity bound |Ret| < r]s^^'^/4:k, (note that ry = 1 in 
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FIG. 6: Using tiie UFD set as input, we sliow tlie decompo- 
sition of Roy equations into the subtracting term ST, kernel 
term KT, and driving term DT for tiie SO, P and S2 waves. 
Note the different scales used on each plot. 



FIG. 7: Using the UFD set as input, we show the decomposi- 
tion of GKPY equations into the subtracting term ST, kernel 
term KT, and driving term DT for the SO, P and S2 waves. 
Note the different scales used here and in Fig. (6] 



the elastic region). For comparison, we present in Fig. [7] 
the same decompositions for the GKPY equations. Note 
the very different scales on both sets of Figures. 

As can be seen in Fig.[6l the ST{s) and KT{s) terms in 
Roy equations become huge at higher energies and suffer 
a large cancellation against each other. This cancellation 



is particularly strong for the SO wave, where, for a suffi- 
ciently large energy, both terms are much larger than the 
unitarity bound. For instance, they are larger by roughly 
a factor of four at 750 MeV, and of eight at 1100 MeV. 

In contrast, as seen in Fig. [7]for the GKPY equations, 
Eq. p^ . the ST terms are constant, and in fact much 
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smaller than the KT{s) terms, which are clearly the dom- 
inant ones. Therefore, no big cancellations between any 
two terms are needed in order to reconstruct the total 
real part of the amplitude. Moreover, we have checked 
that the high energy part, which has been parametrized 
by means of Regge theory, corresponds to somewhat less 
than half of the total DT{s) contribution. Therefore, 
although the DT{s) terms in the GKPY equations are 
larger than in Roy equations due to the fact that there 
is one subtraction less, the contribution coming from the 
amplitudes above 1420 MeV is still small compared with 
the dominant term KT(s). Thus, the high energy behav- 
ior is still well under control. 

Note that, to keep the plots clear, we have only pro- 
vided central values for the moment. In the next section 
we will provide the total uncertainties (the uncertainties 
of each separated contribution were presented in a con- 
ference [33] using a very preliminary UFD set). For our 
purposes it is enough to remark that uncertainties follow 
a similar pattern to these central values. In particular, 
the ST(s) term in Roy equations for scalar waves has a 
large uncertainty due to the poor experimental knowl- 
edge of the Qq scattering length, that becomes larger 
and larger, proportionally to s, as the energy grows, be- 
coming dominant above roughly 450 MeV. In contrast, 
since the GKPY ST term is constant and there are no 
large cancellations, the resulting GKPY equations have 
a much smaller uncertainty in that region. Actually, the 
errors for the GKPY equations in the three waves come 
almost completely from the KT{s) terms. At low ener- 
gies, the effect is reversed and Roy equations provide a 
much more stringent constraint than GKPY. Therefore, 
and as we will show next, they become complementary 
ways of checking our data parametrizations at different 
energies. 



F. Consistency check of Unconstrained Fits 

In order to provide a consistency measure for our 
parametrizations with respect to the dispersive relations 
and sum rules presented in the previous sections, we will 
make use (as we did in previous works) of a quantity simi- 
lar to an averaged / {d.o.f .) distribution. In particular, 
we can consider that a dispersion relation i is well satis- 
fied at a point s„ if the difference A^, defined in Eqs. (IT^ . 
([13]), (fT4| and (fT8|) . is smaller than its uncertainty i^A^. 
Thus, when the average discrepancy verifies 



1 



\p / Ai(s„) 

number of points -^-^ \(5Ai(s„) 



< 1, 



(19) 



we consider that the corresponding dispersion relation is 
well satisfied within uncertainties in the energy region 
spanned by the points s„. In practice, the values of slJ"^ 
are taken at intervals of 25 MeV between threshold and 
the maximum energy where we study each dispersion re- 
lation (1420 MeV for FDR and 1115 MeV for Roy and 



GKPY equations). In addition, we have added a point 
below threshold at s = 2M^ for the Fqo and Fo+ FDRs. 

Similarly, we define discrepancies for the sum rules in 
Eqs. (HH) and dUl), as follows: 



(20) 



In order to calculate the uncertainties (SAi(s„), 51 ^ SJ, 
we have followed two approaches: On the one hand we 
have simply added in quadrature the effect of varying 
each parameter independently in our parametrizations 
from Pi to Pi ± Spi. The errors are symmetric since, in 
order to be conservative, we have always taken the largest 
variation as the final error when changing the sign of Spi. 
This is rather simple but does not take their correla- 
tions into account. On the other hand, we have also es- 
timated the uncertainties using a Monte Carlo Gaussian 
sampling (s^l of all CFD parameters (within 6 standard 
deviations). The uncertainties are then slightly asym- 
metric, corresponding to the independent left and right 
widths of the generated distribution for 10^ events. This 
is, of course, much more time consuming, although in this 
way we can keep part of the correlations in the results. 
However, we have checked that both methods yield very 
similar results, because the errors coming from each indi- 
vidual parameter are small and the number of parameters 
is large. The difference between using one method or an- 
other is almost negligible [11] and thus, for simplicity, we 
are providing numbers and figures with the first one, that 
would be much easier to reproduce should someone use 
our parametrizations. 

In Table U we show the averaged squared discrepan- 
cies 3f that result when we use the UFD set described in 
Sects. inland IIIII We are showing these discrepancies up 
to two different energy regions, 932 MeV and 1420 MeV 
for FDRs, and up to 992 MeV and 11115 MeV for both 
Roy and GKPY equations (note that we have kept the 
same definition of energy regions as in KPY08, so that 
we can compare easily with the results obtained there). 
Let us remark that these discrepancies are "squared dis- 
tances" , similar to a x^, and so we will abuse the language 
and talk about average "standard deviations" , which cor- 
respond to the square root of . Still, one has to keep in 
mind that these dispersion relations have not been fitted 
yet. 

Let us first concentrate in the low energy part below 
932 MeV or 992 MeV. We can observe that FDRs are 
reasonably well satisfied: discrepancies are never beyond 
1.3 standard deviations. Roy equations are also well sat- 
isfied, with a discrepancy below 1.2 standard deviations. 
However, the GKPY equations arc much more demand- 
ing: The UFD set satisfies the S2 wave equation fairly 
well, but it does not satisfy the SO and P wave relations 
so well. Still, no dispersion relation lies beyond 1.6 stan- 
dard deviations. This is not too bad, given the fact that 
we have not fitted the dispersion relations, but there is 
clear room for improvement. Let us recall that this is just 
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new UFD 


old UFD 


new UFD 


old UFD 


FDRs 


8^/'^ < 932 MeV 


s^''^ < 1420 MeV 


7r TT 


0.31 


0.12 


2.13 


0.29 


77 TT 


1.03 


0.84 


1.11 


0.86 


7* 1 


1.62 


0.66 


2.69 


1.87 


xioy Jijqs. 


s^/^ < 992 MeV 


s^/^ < 1100 MeV 


On 


0.64 


0.54 


0.56 


0.47 


Go 


1.35 


1.63 


1.37 


1.68 


P 


0.79 


0.74 


0.69 


0.65 




s^^^ < 992 MeV 


s^''2 < 1100 MeV 


so 


1.78 


5.0 


2.42 


8.6 


S2 


1.19 


0.49 


1.14 


0.58 


p 


2.44 


3.1 


2.13 


2.7 


Average 


1.24 


1.46 


1.58 


1.97 



TABLE I: Average discrepancies of the unconstrained 
data fits (UFD set) for each dispersion relation. We com- 
pare the results of the parametrization obtained in this work 
(new UFD) with those in KPY08 (old UFD set). The huge 
discrepancies seen in KPY08 for GKPY equations all come 
from energies above ~ 500 MeV . This is the main reason to 
improve our unconstrained SO fit as explained in Sect. IIII CI 

how experimental data satisfy these constraints, there is 
no theory on the UFD set. 

If we now also include the region above 932 MeV for 
FDRs or above 992 MeV for Roy and GKPY equations, 
we find that the agreement deteriorates considerably: 
four relations lie between 1.4 and 1.65 average standard 
deviations, but not beyond that. Fortunately we will get 
much better fulfillment of dispersion relations in all re- 
gions by allowing for a small variation of the parameters 
in the constrained fits to be discussed below. 

Let us also remark that the two sum rules, Eqs. (fT6|) 
and (fT7|) . are satisfied within 1.9 and 0.3 standard devi- 
ations. Even for the first one, this is still a fair agree- 
ment, because, in practice, both of them correspond to 
a one-order of magnitude cancellation between the low 
and high energy contributions to the sum rules, which, 
in these UFD set are determined from uncorrelated data 
fits. 

Also in Table U we show the average discrepancies for 
the old UFD set in KPY08. With regard to FDRs and 
Roy equations, it is evident that the new UFD fit is doing 
worse than the one in KPY08. Nevertheless, one should 
keep in mind that the new SO wave has reduced its un- 
certainty at low energies by somewhat more than 10%, 
because the published NA48/2 data are more precise and 
also because we are discarding the controversial K ^ 2tt 
datum. For that reason, one would have expected the 
averaged squared discrepancies to look now bigger by as 
much as 20 or 30% whenever the SO wave contributes 
significantly to the dispersion relation. With this cor- 
rection in mind, the deterioration is not so significant. 
Nevertheless, we want to insist that this is basically due 
to the new results of NA48/2 and our getting rid of the 
K — > 2tt datum. The data have changed. 

Why do we then claim to have improved the SO wave 
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FIG. 8: Results for Forward Dispersion Relations. Dashed 
lines: real part, evaluated directly with the UFD parametriza- 
tions. Continuous lines: the result of the dispersive integrals. 
The dark bands cover the uncertainties in the difference be- 
tween both. From top to bottom: (a) the ■K%° FDR, (b) the 
■K°-K+ FDR, (c) the FDR for 7* = 1 scattering. The dotted 
vertial line stands at the KK threshold. 

in this work? The answer comes from GKPY equations, 
which, as we already explained, are much more precise 
than Roy equations above roughly 450 MeV for the SO 
wave, given the present experimental input. It is clear 
that the KPY08 UFD parametrization satisfies the SO 
GKPY equation very poorly at any energy and is not sat- 
isfying the low energy P GKPY equation very well. For 
that reason, we have improved the matching and the data 
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selection, so that our new UFD parametrization, that 
will be our starting point for the constrained fits, satis- 
fies GKPY equations much better without spoiling FDR 
and Roy equations. The improvement due to the new 
unconstrained SO wave fit is obvious from Table HI par- 
ticularly in the SO GKPY equation. Up to 1100 MeV, the 
old UFD set from KPY08 had an averaged squared dis- 
crepancy of 8.6, whereas the new UFD set has 2.42. This 
huge improvement on the SO wave has been compensated 
by some deterioration in other relations at high energy, so 
that the averaged discrepancy up to high energies is re- 
duced only from 1.97 to 1.58. Note that the change in the 
inelasticity parameter, that now shows a much bigger dip 
in the 1000 to 1100 MeV region, as shown in Fig. El plays 
a relevant role in this dramatic improvement. This dip 
structure is thus favored by the GKPY equations, some- 
thing that could not be seen with standard Roy equations 
since their uncertainties in that region are huge. We will 
discuss this in detail in Sect. IVII Bl At low energies, the 
average squared discrepancy has been reduced very little, 
from 1.46 down to 1.24. Of course, let us remark once 
again that our uncertainties are now 10-15% smaller in 
the SO wave at low energies, so that the improvement is 
actually bigger than it seems just from the numbers in 
the table. 

Let us mention here that the inclusion of the new terms 
parametrizing a crude dependence on the r] momentum 
above -q-q threshold, help reducing the average squared 
distances by 6%, namely, from 1.68 to 1.58. In particular, 
the averaged squared discrepancies df for the SO GKPY 
equation decrease from 3.02 to 2.42 and for the Foo FDR 
equation from 2.35 to 2.13. 

Up to now, we have studied the overall uncertainties, 
but in Fig. [8] we show to what extent FDRs are satisfied 
by the UFD set, as a function of energy. Of course, the 
best fulfillment is found at lower energies. In Fig. [S] we 
show how the usual, twice subtracted Roy equations are 
satisfied by the UFD set. Here, as we did in Sect. IIVEI 
we denote by "in" what our parametrizations give for 
Re tf^ , whereas we denote by "out" the result of the dis- 
persive representation from Roy equations, namely, the 
subtraction constant terms, plus the kernel terms, plus 
the driving terms in Eq. (jl4p . Finally, in Fig.[TU]we show 
how the new, once subtracted, GKPY equations are sat- 
isfied by the UFD set. We follow the same "in" and "out" 
notation as for Roy equations. 

Comparing Fig. IH] with Fig. [TUJ it is clear that, given 
the present experimental input, the uncertainty band for 
GKPY equations is much smaller than that for Roy equa- 
tions above 450 MeV, whereas the opposite occurs at 
lower energies. Therefore, as we have emphasized repeat- 
edly, the new GKPY equations represent a much stronger 
constraint in the intermediate energy region than stan- 
dard Roy equations. 

In summary, with the new SO unconstrained fit, all 
dispersion relations are satisfied in the different energy 
regions within less than 1.6 standard deviations in the 
low energy regime, and 1.7 including the intermediate 
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FIG. 9: Results for Roy Equations. Dashed lines ("in"): real 
part, evaluated directly with the UFD parametrizations. Con- 
tinuous lines ("out"): the result of the dispersive representa- 
tion. The gray bands cover the uncertainties in the difference 
between both. From top to bottom: (a) SO wave, (b) S2 
wave, (c) P wave. The dotted vertial line stands at the KK 
threshold. 



energies. This is a fairly reasonable fulfillment, given the 
fact that the information about analyticity has not been 
included as a constraint in the UFD description. Never- 
theless, it is obvious that there is room for improvement, 
which is what we will do by obtaining constrained data 
fits in the next section. 
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FIG. 10: Results for GKPY equations. Dashed lines ("in"): 
real part, evaluated directly with the UFD parametrizations. 
Continuous lines ("out"): the result of the dispersive rep- 
resentation. The gray bands cover the uncertainties in the 
difference between both. From top to bottom: (a) SO wave, 
(b) S2 wave, (c) P wave. Note how these uncertainties are 
much smaller above 450 MeV than those from standard Roy 
equations shown in Fig. [9l The dotted vertial line stands at 
the KK threshold. 

V. FITS TO DATA CONSTRAINED BY 
DISPERSION RELATIONS 



In previous works (PY05, KPY08) we had improved 
the consistency of our description of tttt scattering am- 
plitudes by imposing FDR and Roy equations fulfillment 



within uncertainties. As we have just seen in the previous 
section, the GKPY equations provide a much more strin- 
gent constraint in the intermediate energy region than 
standard Roy equations and thus it now makes sense to 
impose the new GKPY equations as an additional con- 
straint in a new set of Constrained Fits to Data (CFD 
set). 



A. Minimization procedure 

Our goal is then to obtain a fit to data, by changing the 
UFD parametrizations slightly, that fulfills each disper- 
sion relation within errors. As we did in Q, we will now 
use the average discrepancies df, defined in Eqs. 
and (j20p . to obtain these constrained fits, by minimizing: 



E 



W^d^ + dj + d'j + 



E 

k 



Pk 



cxp 

Pk 



Spk 



(21) 



where i runs over the three FDRs, the three Roy and the 
three GKPY equations. Here, we denote by p'j^^ all the 
parameters of the UFD parametrizations for each wave 
or Regge trajectory. In this way we force the previous 
data parametrizations to satisfy dispersion relations and 
sum rules within uncertainties. In KPY06 and KPY08 
a common weight of ^ 9 was estimated from the 
typical number of degrees of freedom needed to describe 
the shapes of the output. This value ensured that every 
single dispersion relation was fairly well described by the 
KPY08 constrained data fits up to the matching energy 
used in that work, namely, 932 MeV. 

However we are now considering partial waves up to 
1115 MeV. For most waves, this extension does not alter 
significantly their shape and Wi = 3 is still a good weight. 
Nevertheless, we have less points in the region above 932 
MeV and if we want the fit to give not just a good aver- 
age df, but also a good description for each wave, some 
of these waves need further weight on the high energy 
region, in particular if their UFD df was larger than 2. 
For this purpose, we have increased Wi up to 3.5 for the 
high energy parts of the Fqo, F*^-^*^^', as well as 4.2 for 
the GKPY P wave in the whole energy region. Finally, 
we have increased Wi up to 7 for the high energy part of 
the SO GKPY equation. The latter was to be expected, 
since in this region there is a lot more of structure, both 
in the phase and inelasticity, due to the presence of the 
/o(980). These values are not arbitrary, since they have 
been obtained by increasing each Wi gradually, start- 
ing from 3, until the df are below or very close to one 
uniformly throughout the whole energy range, for all dis- 
persion relations obtained from the constrained fit. This 
uniformity is very relevant to avoid dispersive constraints 
being badly satisfied in some small energy region despite 
the averaged df still remaining below 1. 

Before proceeding further, let us recall that, strictly 
speaking, the quantity that we minimize in Eq. (1211) is 
not a x^, but that each individual df is a measure of how 
well each dispersion relation is satisfied. 
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Roy Eqs. 


GKPY Eqs. 


Roy Eqs. 


GKPY Eqs. 


with UFD 


with UFD 


with CFD 


with CFD 


112 ±24 


120 ± 30 


83 ±32 


85 ±34 


\/s|2 189 ±11 


200 ±6 


200 ± 10 


201 ± 5 



TABLE II: Adler zero positions y/FX, in MeV, for the SO and 
S2 waves, obtained from Roy or GKPY equations using either 
the parametrizations form the UFD or CFD sets. 



B. Variation of the S2 Adler zero 

As we have seen in Section IIII C[ in the parametriza- 
tion of each scalar wave we exphcitly factorized a zero 
in the subthreshold region. These are the Adler ze- 
ros required by chiral symmetry constraints (30| . Ac- 
tually, we fixed them to ee ^ A/2/2 ~ 99 MeV and 

y^s^ EE y/2M^ ~ 197 MeV, which are their current al- 
gebra values (leading order ChPT). Of course, once these 
UFD parametrizations are used inside the SO and S2 Roy 
or GKPY equations we can also obtain the dispersive re- 
sult for the SO and S2 Adler zeros, that we provide in 
Table ini 

In order to determine the positions of Adler zeros bet- 
ter when making constrained fits in KPY08, we allowed 
them to change within the dispersive uncertainties ob- 
tained from the UFD set. However, in this work we will 
not insist on Zq/V^ reproducing the SO wave Adler zero 
very precisely. The reason is that, as we see in Table HIl 

the uncertainties in ^ s^ obtained either from Roy or 

GKPY equations are huge, and setting zq free introduces 
a spurious and extremely correlated source of error. In 
addition, in KPY08, the zq central value moved in the 
wrong direction [4l']. In addition, as already explained 
in Sect. IIII G[ the SO wave Adler zero lies close to the 
border of the conformal circle, i.e., w{sj^) ~ —1, where 
the conformal expansion coverges very slowly. We simply 
have to accept that our SO wave conformal expansion is 
not very accurate around the Adler zero. Of course, this 
is irrelevant for the integrals in the physical region and 
has a negligible infiuence in the set of constrained fits we 
will obtain next. 

In contrast, the S2 Adler zero obtained from the dis- 
persive representation moves very little from its current 
algebra value and its uncertainty is rather small. The 
reason for this difference in uncertainties is, for a good 
part, that the SO wave Adler zero lies very close to the 
left cut, whereas the S2 Adler zero is not so far from 
threshold and is quite well determined when data is used 
as input of either Roy or, even better, GKPY equations. 
For that reason, we still allow the S2 Adler zero to vary 
when making the constrained fits, using as a starting 
point the weighted average of the values obtained from 
the UFD set inside Roy and GKPY equations, namely. 




= 197.7 ±5.1 MeV. 
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< 992 MeV 
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so 




0.23 




0.24 


S2 




0.12 




0.11 


P 




0.68 




0.60 


Average 


0.22 


0.28 



TABLE III: Average discrepancies of the Constrained Fits 
to Data (CFD) for each dispersion relation. 



C. The Constrained Fits to Data (CFD) 

The resulting parameters for the CFD are gathered in 
the tables of the Appendix A. It is reassuring to observe 
that, except for the SO wave at intermediate energies, 
the values of the parameters do not change much from 
the UFD to the CFD sets, as could be expected, since, 
as we saw in Table HI the UFD fulfillment of dispersive 
constraints only needed some improvement, but not a 
radical change. In particular, the GKPY equation for 
the SO wave is very well satisfied in the CFD at the ex- 
pense of an average change of 0.82 standard deviations 
in the high energy parameters and almost no change in 
the low energy ones. Certainly, most of this change is 
concentrated in the parameters c and ei in Eqs. (|10p and 
dm . We will discuss below that the resulting phase after 
this change still describes the phase shift and inelasticity 
data fairly well, but tends to make the /o(980) somewhat 
wider. The D2 wave is the one that deviates most from its 
unconstrained parametrization, but its parameters are, 
on average, within 1.4 standard deviations of their UFD 
value. This could be expected, as was already commented 
in our previous works [H, since, together with the SO 
at high energy, it is probably the one where data has 
the worst quality. The parameters of the other waves, or 
those of the Regge parametrizations, do not deviate — on 
the average — beyond 0.6 standard deviations from their 
UFD values. In Table IXU in Appendix ID] we provide the 
SO, P and S2 phase-shifts that result from using the CFD 
set inside the dispersive representation. 

In Table IIIII we list the averaged discrepancies that re- 
sult when we use the constrained fits (CFD) inside the 
dispersion relations. Let us remark that all discrepan- 
cies are now below one, and very similar both for the low 
energy region and also when including the high energy re- 
gion. This shows a remarkable average consistency and 
homogeneity for this new set of data parametrizations. 
Let us recall that we only constrain our fits to satisfy 
dispersion relations up to 1420 MeV for FDR and 1115 
MeV for Roy and GKPY equations. Consequently, we ex- 
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FIG. 11: Results for forward dispersion relations. Dashed 
lines: real part, evaluated directly with the CFD parametriza- 
tions. Continuous lines: the result of the dispersive integrals. 
The dark bands cover the uncertainties in the difference be- 
tween both. From top to bottom: (a) the tt^tt" FDR, (b) the 
7rV+ FDR, (c) the FDR for It = 1 scattering. 

pect the dispersive representation to be somewhat worse 
satisfied in the region near the maximum energy under 
consideration. This is indeed observed since the average 
squared discrepancies are somewhat smaller below 1 GeV 
than up to the maximum energy, where we usually find 
the point satisfying the dispersion relations worse. 

Furthermore, as already commented, the updated se- 
lection and treatment of the SO wave data has decreased 



the SO wave uncertainties by roughly 10 to 15%. This 
means that the consistency shown by the average dis- 
crepancies in Table IIIII is even better than it looks when 
comparing with similar results given in KPY08 for FDR 
and Roy equations, since we are getting a very good con- 
sistency with slightly smaller uncertainties. 

As we did for the UFD set, we now show in Figs. [11] 
[m and [13 how well the CFD set satisfies FDR, Roy and 
GKPY equations respectively. The improvement in the 
consistency of the CFD set over the UFD is evident by 
comparing these plots with their UFD counterparts in 
Figs. [51[n| and [ini 

Finally, the two sum rules in Eqs. (flSl) and p7|) are also 
remarkably well satisfied, within 0.93 and 0.1 standard 
deviations, respectively. In particular, the 1.9 standard 
deviations for the sum rule in Eq. PT|) using the UFD 
set are reduced dramatically, and this implies now a two 
orders of magnitude cancellation between the low and 
high energy contributions. 



VI. THRESHOLD PARAMETERS AND ADLER 
ZEROS 



Apart from the additional GKPY equations, the main 
novelty of this work is the SO wave improvement, both in 
its parametrization and data analysis. Thus, naively, one 
may not expect a big variation in the low energy part of 
the other waves with respect to previous works. 

However, let us recall that, as we did in KPY08, we cal- 
culate most threshold parameters from sum rules. Thus, 
the changes in the SO wave can also affect the calculation 
of these low energy parameters for other waves. In par- 
ticular, when using sum rules with one subtraction, the 
intermediate energy part of our parametrizations, now 
constrained by GKPY equations, also plays a relevant 
role in our final results. In this section we will thus recal- 
culate all these threshold parameters with the new CFD 
set. Actually, we will find that not only the SO wave, but 
the D wave threshold parameters suffer sizable modifica- 
tions. 

Finally, in previous works we did not use the disper- 
sive or sum rule techniques to determine with precision 
the position of Adler zeros, which are required by chiral 
symmetry in the subthreshold region of the SO and S2 
waves, and are therefore of interest for Chiral Perturba- 
tion Theory. Also in this section we will determine them 
using the Roy and GKPY equations with the CFD set as 
input for the integrals. 



A. Sum rules for threshold parameters 

We list in Table IIVI the values of threshold param- 
eters for all the partial waves we considered in this 
analysis, namely: SO, S2, P, DO, D2 and F. In addi- 



tion, we provide values for Oq 
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FIG. 12; Results for Roy equations. Dashed lines ("in"): 
real part, evaluated directly with the CFD parametrizations. 
Continuous lines ("out"): the result of the dispersive rep- 
resentation. The gray bands cover the uncertainties in the 
difference between both. From top to bottom: (a) SO wave, 
(b) S2 wave, (c) P wave. 



6q^\m]^) — (5q^''(M|-), since these parameters are of rel- 



evance for pion atoms, scalar threshold parameters, and 
kaonic decays. In the second and third columns, we pro- 
vide the results from the UFD and CFD sets. We al- 
ready commented that the CFD parametrizations change 
only very slightly compared to the UFD, and this is well 
corroborated by the fact that all the UFD and CFD re- 
sults on Table HV] are compatible with one another within 
roughly one standard deviation. 



s (MeV) 

FIG. 13: Results for GKPY equations. Dashed hues ("in"): 
real part, evaluated directly with the CFD parametrizations. 
Continuous lines ("out"): the result of the dispersive rep- 
resentation. The gray bands cover the uncertainties in the 
difference between both. From top to bottom: (a) SO wave, 
(b) S2 wave, (c) P wave. Note how these uncertainties are 
much smaller above 450 MeV than those from standard Roy 
Eqs. shown in Fig. [12] 

In the fourth column, we use the very reliable CFD set 
inside several sum rules, that we detail next only very 
briefly, since they had already been given in detail in 
KPY08. First, we use the well known Olsson sum rule: 



Imf(^'-i)(5,0) 
4M5 .(.-4M2) ' ^''^ 
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which is dominated at high energies by the p-Regge ex- 
change, and can thus have only one subtraction. Apart 
from the normalization, this is just the FDR in Eq. ([T3l) . 
but evaluated at threshold. 

Next, for £ > 1, we use the Froissart-Gribov represen- 
tation: 



at 



b,= 



V^r(^ + i) 

4M„r(£ + 3/2)74^2 



2Af„r(£ + 3/2)74^2 "I ( 
(^-K l)ImF(s,4Af2) 



4M2). 



(23) 



with IxnF^^^g = {d / d cos 9 s)lm F , where cos 6*3 is the an- 
gle between the initial and final pions. For amplitudes 
with fixed isospin in the t channel, an extra factor of 2 
(due to identity of particles) has to be added to the left 
hand side of the equation above. 

In addition, we use the following sum rule that we de- 
rived in [l|: 



bi 




1 



lmi^(-^*=°'(s,0) 



.(S-4M2) 



Imi^(^t=i)(s^O) 
ImF(^*=2)(s,0) 



(24) 



together with another two sum rules, derived in in- 
volving either the SO and S2 slopes: 



6^"' + 26^'' 



lim P.R/ d 



(25) 



, 6Af^(2s'-4M2)ImFoo(s') 



4M^+ ■ JiM^ s'is' + s ~ 4M2)(s' - 4Af2)(s' _ s) ' 

or the S2 slope parameter and the P wave scattering 
length: 



3aW+6f = 



lim P.P./ ds 



, 4A//,(2s'-4A/2)ImFo+(s') 



(26) 



,^4Af5+ ■ 'J^Ml s'(s' + s-4Af2)(s'-4M2)(s'-s)' 

Note that, as explained in Q, the limits above are to 
be taken for s > AM^. In practice, for the value of oi 
we simply use its Froissart-Gribov representation and we 
are left with a sum rule representation for both b'^'^ and 

^(2) 

The results for all these sum rules are listed in the 
fourth column of Table IIVI 



The fifth column, that contains what we consider our 

best values, is obtained as follows: For 2a'i^^ — 5a'i^\ 6q°\ 
(2) 

&Q , ai and bi, we take the average between the sum 
rules above and the direct value of the CFD set, since 



UFD CFD Sum rules with CFD Best values KPY08 values 

0.218 ± 0.009 0.221 ± 0.009 0.220± 0.008" 0.223 ± 0.009 

4^' -0.052 ± 0.010 -0.043 ± 0.008 -0.042± 0.004'= -0.044± 0.004 



% - % 0.270 ± 0.009 0.264 ± 0.009 0.262 ± 0.006'= 0.267 ± 0.009 

24"' - 54^' 0.696 ± 0.054 0.657 ± 0.043 0.648± 0.016" 0.650 ± 0.015 0.668 ± 0.017 



,(0) 
) "0 

ly — 5aQ 

S^°\M]f) - S^iMjc) 47.4 ± 0.9° 47.3 ± 0.9° 47.3 ± 0.9° 50.9 ± 1.2° 

b^°' 0.276 ± 0.007 0.278 ± 0.007 0.278± 0.008'' 0.278 ± 0.005 0.290 ± 0.006 

b^g^ -0.085 ± 0.010 -0.080 ± 0.009 -0.082± O.OOl"* -0.082 ± 0.004 -0.081 ± 0.003 

ai(xlO^) 37.3 ± 1.2 38.5 ± 1.2 37.7 ± 1.3*" 38.1 ± 0.9 38.1 ± 0.9 

6i(xlO^) 5.18 ± 0.23 5.07 ± 0.26 6.0 ± 0.9\ 5.48 ± 0.17'= 5.37 ± 0.14 5.12 ± 0.15 

4°'(xl0'') 18.7 ±0.4 18.8 ± 0.4 17.8 ± 0.3*" 17.8 ± 0.3 18.33 ± 0.36 

4^'(xl0'') 2.5 ± 1.1 2.8 ± 1.0 1.85 ± 0.18** 1.85 ± 0.18 2.46 ± 0.25 

6^°'(xlO'') -4.2 ± 0.3 -4.2 ± 0.3 -3.5 ± 0.2'' -3.5 ± 0.2 -3.82 ± 0.25 

4^'(xl0'') -2.7 ± 1.0 -2.8 ± 0.8 -3.3 ± O.l'' -3.3 ± 0.1 -3.59 ± 0.18 

a3(xl0^) 5.2 ± 1.3 5.1 ± 1.3 5.65 ± 0.23' 5.65 ± 0.21 6.05 ± 0.29 

63(xl0-'^) -4.7 ± 2.6 -4.6 ± 2.5 -4.06 ± 0.27'' -4.06 ± 0.27 -4.41 ± 0.36 

TABLE IV: Threshold parameters in the customary M-n- = 1 units and the 5q''''(M|-) — 5q^''(A4"|-) phase difference. The 
values in the second and third columns are obtained directly from the UFD and CFD parametrizations, respectively. 
The fourth column is obtained using the CFD set inside sum rules: '^from Eq. (f^ . ''from Eq. (^5)1 . ''from Eq. (f^. 
''from Eqs. (|25|) and ((26}. In addition, for the scalar scattering lengths'^ best values, we have re-fitted their CFD 
values constrained to satisfy the Olsson sum rule, Eq. ([22]), which is also used to obtain the best value for their 
difference and its uncertainty, Eqs. (|27p and 
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they are basically independent. However, for the DO, 
D2 and F waves, in order to stabilize the fits, we had 
already constrained the value of the threshold parameters 
by means of the Froissart-Gribov representation in the 
UFD set (see Hence, in those cases, it makes no 

sense to average either the UFD or CFD direct result with 
the Froissart-Gribov for 0^2'^ , 0.^2^, b^2^ and 03, which is 
therefore considered our best result. The only exception 

are 62 ^^'^ ^3: since those values were not constrained 
in the initial UFD, but their uncertainty from the CFD 
is an order of magnitude larger than from the sum rule, 
which value is the one we quote as the best one. 

Let us remark that the SO and S2 scattering lengths, 
which are of special interest for ChPT, are refined by 
re-fitting them again to the CFD direct results and the 
Olsson sum rule simultaneously. Obviously, the resulting 
errors are strongly correlated and the corresponding cor- 
relation ellipse is shown in Fig. [141 The uncertainties can 
be uncorrelated by using two new variables, x, y defined 
as: 

4°' === 0.220 + 0.130 X + 0.337 y, 
af'^ = -0.042 -0.337 a; + 0.130?/, 
4°^ - = 0.262 + 0.467 a; + 0.206 y, 
a;=0±0.076, y = ± 0.023, (27) 

which give the numbers listed in the tables as our "Best 
values" : 

= 0.220 ± 0.008, (28) 

= -0.042 ± 0.004, 

- = 0.262 ± 0.006, 

in units of Mjr . 




0.205 0.210 0.215 0.220 0.225 0.230 0.235 



FIG. 14: The one and two standard deviation ellipses (thick 
and dashed lines, respectively) in the {a^'^ , a''Q ^ ) plane. The 
square covers the uncertainties of our best results in Eq. (|28|l , 
obtained from the uncorrelated expressions in Ea. (|27|) . 

For the sake of comparison, we list in the sixth column 
our results from KPY08, where we did not impose the 
GKPY equations nor did we add the several improve- 
ments to the amplitudes and the data implemented in 
this work. Note that the low energy parameters are quite 
consistent with our previous results, i.e. many central 



values lie within one standard deviation of our KPY08 
results, and most of them overlap within one standard 
deviation. There are, of course, the expected exceptions: 
first, the 6^°^ (M|) - S^^'' (M|) central value changes by 3 
standard deviations, mostly due to the fact that we have 
discarded here the controversial K — tttt datum. Next, 
the SO slope ^q"^ changes by two standard deviations, and 
this is mostly due to the inclusion of the isospin viola- 
tion correction in the low energy K^^ data. One could 
have expected that the scattering length a'^^ may have 
suffered a large shift for the same reason, but it has only 
decreased about a third of a standard deviation. Hence, 
most of the change due to the isospin correction is 
concentrated on the slope parameter. In addition, as we 
already anticipated, both D wave scattering lengths have 
decreased by roughly 2 standard deviations. 

Although it will be commented in detail in the dis- 
cussion section, let us note that these new results are in 
much better agreement with the results in [ill] , than were 
those in KPY08. 

As commented in Sec. |Vl we can also check here that 
the new uncertainties are slightly smaller, but only by 
10-15%, than in KPY08, due to: discarding the K —i' tttt 
conflicting input, keeping the SO Adler zero fixed, and the 

more precise NA48/2 published data. The Cq"' — Cq^'' un- 
certainty in has decreased by almost 50%, although 
this is not only due to our improvement of the SO wave, 
but mainly to the fact that we are now calculating it 
differently, using Eqs. (P?]) . 



B. Determination of Adler zeros 

As already explained, chiral symmetry requires the ex- 
istence of zeros in the amplitude close to s = for the 
scalar waves SO and S2 . We have explicitly factorized 
them in our amplitudes at = Zq/2 and = 2z|, see 
Eqs. © or (jAll) and (|A5|) . As a starting point, we have 
first fixed them to the ChPT leading order estimate by 
setting zq — Z2 — Mt^ for the UFD parametrizations. We 
then used these parametrizations inside Roy or GKPY 
equations to recalculate the position of these Adler ze- 
ros, which were listed in the first two columns of Table 

In previous works, we allowed the zo and Z2 parame- 
ters to change in the CFD set, expecting them to be ac- 
curately fixed by imposing the dispersion relations. Un- 
fortunately, as discussed in Sect. IV 51 this does not work 
for the SO wave. The reason is that its Adler zero is 
very close to the left cut, in a region where, on the one 
hand, neither Roy nor GKPY equations provide a precise 
determination of the zero position (see Table HII and, on 
the other hand, the conformal expansion converges badly. 
For that reason, we have simply kept the SO parameter zg 
fixed to Airr both on the UFD and CFD sets. Being so far 
from the threshold region, this effect is irrelevant inside 
the dispersive integrals. Thus, only the S2 Adler zero 
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is allowed to change when obtaining the CFD set, but 
only within the UFD uncertainties obtained from Roy 
and GKPY equations. 

In this section we go one step further and we finally 
provide, in the last two columns of Table HIl the value 
of the SO and S2 wave Adler zeros obtained when the 
CFD set is used inside Roy and GKPY equations. The 
CFD SO zero is closer to its expected position (around 80 
MeV) than the UFD result, but note that the uncertainty 
gets worse because of this displacement towards the left 
cut. In summary, we do not have enough precision to pin 
down the location of this SO Adler zero accurately. 

In contrast, the S2 Adler zero is determined quite pre- 
cisely by GKPY equations (and to a lesser extent by Roy 
equations), and the resulting zi parameter, if allowed to 
vary, is almost identical to its UFD determination. Thus, 
as explained in Section FyBl we have allowed zi\/2Xo vary 
within the weighted average between the GKPY and Roy 
equation results of the UFD set. The resulting Adler 
zero, when read directly from the CFD parametrization 

is \J^^ — Z2\/2 — 201 ± 5 MeV, which is almost iden- 
tical to the values obtained by using the CFD set inside 
Roy or GKPY equations — listed in Table |lll This con- 
firms that it is correct to identify the Adler zero with the 
Z2-\/2 term in our S2 wave conformal parametrization. 



VII. DISCUSSION 

First of all, we want to remark that ours is just a data 
analysis, and we are not predicting the value of any ob- 
servable, just determining them from experiment. In con- 
trast to other approaches [ll|, we are not solving FDR, 
Roy or GKPY equations, but just imposing them as con- 
straints on the data analysis. Actually, all these equa- 
tions have been obtained with several approximations, 
for instance, they are obtained in the isospin limit, and 
we only expect them to describe the real world up to 
some uncertainty of the order of 3%. In addition, all Roy 
equations studies we are aware of — including this one — 
neglect any inelasticity to four or more pion states below 
the two kaon threshold. This is certainly a very small 
effect, but is nevertheless an approximation. 

Being a data analysis, our parametrizations change 
when the data changes. In particular, in this work we 
have updated the NA48/2 data Q with their final results 
[20j . which have smaller uncertainties. In addition, we 
have incorporated the threshold-enhanced isospin correc- 
tion in 22] to all K^^ data. Moreover, we have discarded 
the controversial K ^ tttt datum [2^. Furthermore, the 
increased precision provided by the once-subtracted dis- 
persion relations that we have introduced in this work, 
requires an improved parametrization with a continuous 
derivative matching. This additional constraint and the 
requirement that the output of the dispersion relations 
should satisfy the elastic unitarity bound — which is auto- 
matic in the input parametrizations — , has made us also 



add an additional parameter to the SO wave parametriza- 
tion at low energies. As we will see below, the SO wave 
parametrization at intermediate energies favors the "dip- 
scenario" for the SO inelasticity between 1000 and 1100 
MeV. In this discussion section we will show in detail 
the new CFD set, particularly the SO wave, comparing it 
to other works, and we will discuss the consequences of 
these modifications. 



A. The new CFD SO wave 

In Fig. [15] we show the resulting CFD SO wave from 
threshold up to 1420 MeV, versus the data from different 
sets in the literature [13, HI] • Note the smooth matching 
at 850 MeV and the kink at KK threshold. This is in 
contrast with our old KPY08 results, already shown in 
Fig. [31 which have a spurious kink at the matching point 
(932 MeV in that work), and a much less pronounced 
kink at KK threshold. The difference between the UFD 
and CFD SO wave phase shift at low energies, which we 
showed in Fig. [21 is almost imperceptible. 

To ease the comparison of this CFD results with the 
UFD set for all energies, we have plotted their central 
values together in Fig. [THl It can be noted that the 
change above KK threshold is again almost impercep- 
tible up to 1200 MeV. The only sizable differences be- 
tween the phase of the UFD and CFD parametrizations 
are above 1200 MeV, where our parametrizations are less 
reliable since Roy and GKPY Eqs. only extend up to 
1115 MeV, and on the sharp phase rise in the 900 MeV 
to 2mK = 992 MeV region due to the /o(980) resonance, 
which is clearly less steep in the CFD case than in the 
UFD. The latter is one of the reasons why the CFD solu- 
tion satisfies GKPY equations well within uncertainties, 
but the UFD lies somewhere around 2 standard devia- 
tions away (see Tables |T] and IIIIl respectively) . 

In addition, we are also showing in Fig. [THl the results 
from 11], which are in good agreement with ours, but lie 
slightly lower only above 550 MeV (see discussion below) . 
Actually, our CFD solution does not show the "hunch- 
back" between 500 and 900 MeV seen in KPY08, as al- 
ready shown in Fig. [3l 

Concerning the SO inelasticity, we show in Fig. [T71 the 
difference between the UFD and CFD sets. It can be no- 
ticed that the difference lies essentially within the uncer- 
tainties (gray area), although the "dip" structure above 
1000 MeV becomes even deeper in the CFD set. Finally, 
in Fig. [T51 we show the CFD inelasticity versus all the 
existing experimental data. 

Since the UFD set was already providing a good de- 
scription of the inelasticity data obtained from tttt — >■ tttt 
experiments, as shown in Fig. [H so it does the CFD. For 
the same reason, it also fails to reproduce the inelasticity 
data from tttt KK, as we had already shown for the 
UFD case in Fig. [SI Note that this is due to the fact that 
both our UFD and CFD solutions show a "dip" struc- 
ture between 1 and 1.1 GeV, which is seen in the data 
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FIG. 15: The new constrained fit (CFD) for the SO wave versus the existing phase shift data from [27l. [2^. The dark band 
covers the uncertainties, 




FIG. 16: Comparison between the phase of the constrained (CFD) and unconstrained (UFD) Fit to Data for the SO wave. We 
also plot the phase from the Roy equations analysis in [Tlj . 

coming from tttt — > tttt, but not in those coming from allowing for a large deviation from the phase shift data. 
TTTT — KK. This is a longstanding problem (see [2l| and 
references therein) that we will address in the next sub- 
section, showing that the "non-dip" scenario is not able 
to satisfy the dispersive representation well even when 
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FIG. 17: Comparison between the UFD and CFD SO wave 
inelasticity. The gray area corresponds to the CFD uncer- 
tainty. A similar size area should be associated to the UFD 
result, but for clarity we only show its central value. Note the 
dip structure between 1 and 1.1 GeV. 




FIG. 18: CFD SO wave inelasticity versus experimental data. 



B. SO inelasticity: the "non-dip" scenario is 
disfavored 



In order to show how much the "non-dip" scenario is 
disfavored, we will first repeat the same procedure of this 
whole paper, but starting from the SO inelasticity fitted 
to the "non-dip" data, as shown in Fig. [191 while keeping 
the same UFD parametrization for all other waves and for 
the SO phase. We will refer to this set as "ndUFD" . The 
resulting averaged discrepancies 3f are relatively similar 
to those in Table U for our UFD, except for the SO wave 
GKPY equations up to ^/s < 1100 MeV, whose averaged 
df rises from 2.42 to 4.77. This already disfavors the 
"non-dip" scenario. 

Of course, the "dip scenario" UFD set was not doing 



FIG. 19: SO wave inelasticity versus the "non-dip" tttt — ^ KK 
data. We first show the "ndUFD" set obtained from a fit to 
these "non-dip" data. Next, we show the "ndCFD" set ob- 
tained with enlarged errors to try to fulfill dispersion relations. 
This constrained fit satisfies the dispersive constraints better 
but does not describe these "non-dip" data, coming closer to 
the best CFD set, which actually describes the alternative 
"dip" data from tttt — >■ tttt. 



very well either, but we were able to improve it by con- 
straining the fit to data with dispersion relations, i.e., the 
CFD set. One could wonder if a similar quality fit can 
also be obtained by imposing the dispersive constraints, 
but starting from the "ndUFD" . Thus, we followed again 
the procedure described in previous sections, but to ar- 
rive now to a "ndCFD" set. Surprisingly, the SO inelas- 
ticity barely changes, but the improvement comes from a 
bigger variation of the phase in the two-kaon subthresh- 
old region. The resulting average discrepancies c?f come 
in general larger than for our CFD set, sometimes by a 
factor of two, but still below 1. This may look like an 
agreement, but one should not be misguided now by these 
relatively low averaged (If because, contrary to the CFD 
set where discrepancies are below 1 uniformly over the 
whole energy region, for the ndCFD set they are larger 
in the /o(980) resonance region. 

In particular, in the interval between 950 and 1050 
MeV, for the CFD set, the GKPY S0_ equation have 
cP = 1.02, whereas the ndCFD set has (f = 3.49. This 
averaged discrepancy is unacceptable now, since this time 
we are using the dispersion relations as constraints of our 
fits. In addition, the crossing sum rule in Eq. (I16p grows 
to dj = 2.0. 

Furthermore, as we show in Fig. 1201 in the region from 
900 MeV up to KK threshold, the resulting phase of 
this ndCFD scenario lies above all data points with a 
X^/#Points — 3.4, which is a very bad fit. In contrast, 
the CFD set has x^/#points = 0.98 in this region and 
is just a small modification from the UFD phase, which 
has x^/#points = 0.63. Moreover, the ndCFD parame- 
ters lie far from the original ndUFD ones, with the c pa- 
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FIG. 20: Comparison of the UFD, CFD and ndCFD solu- 
tions for the SO phase in the 850 to 1050 MeV region. Note 
that the ndCFD parametrization is largely inconsistent with 
data, despite the fact we are plotting the PY05 averaged data, 
that includes our estimations of the large dominant system- 
atic uncertainties. 



rameter more than 6 standard deviations away from its 
ndUFD value. These numbers clearly show the incom- 
patibility of the ndCFD set with the SO wave tttt — > tttt 
phase-shift scattering data. This disagreement cannot be 
mended by adding systematic uncertainties, since in this 
region we had already included large systematic uncer- 
tainties (see KPY08 and PY05 for details) and all points 
have total uncertainties of more than 10 degrees. 

One could wonder if our minimization procedure, that 
was good enough to reach d? < 1 for the dip scenario, 
is badly tuned for the "non-dip" one. This, of course is 
the role of the Wi weights in Eq. ([?!]) . For this reason 
we have repeated the above procedure adding additional 
weight to the GKPY SO wave equation above 900 MeV. 
The resulting "ndCFD2" yields = 2.06 for the GKPY 
SO equation. Besides, the crossing sum rule in Eq. 
is also df = 1.43. Although they still disfavor this solu- 
tion, these numbers by themselves are not too bad. How- 
ever, the phase shift data between 950 and 1050 MeV has 
X^/#points — 5.9, so that it is described even worse than 
with the previous ndCFD. 

Since we cannot fix the dispersive constraints without 
spoiling the data phase description, we have allowed, as 
a final check, for larger errors in the inelasticity parame- 
ters of the "non-dip" scenario, and applied the dispersive 
constraints. In so doing, we can obtain dj < 1 uniformly 
over all energy regions for all GKPY equations except for 
the SO wave between 950 to 1050 MeV, for which we ob- 
tain 1.42. However, the central value of the inelasticity 
for the resulting constrained "non-dip" fit starts devel- 
oping a dip as seen in Fig. [121 Therefore, we describe 
neither the "non-dip" nor the "dip" scenario. 

In conclusion, the non-dip scenario, even when con- 
strained with dispersion relations, is not able to describe 



the data and simultaneously satisfy Forward Dispersion 
Relations, Roy and GKPY equations, plus certain cross- 
ing sum rules. 



C. Comparison with other works. 



The results listed in Table ITVl for threshold parameters 



are remarkably compatible with the predictions of 
using Chiral Perturbation Theory and Roy equations 



Hi 



.(0) 



0.220 ±0.005, 



,(2) 



-0.0444 ±0.0010. 



The agreement with that reference has also improved a 
great deal since the d^°\M'^) - S^iM]^) = (47.3 ±0.9)° 
value, obtained directly from our CFD set, is now com- 
pletely consistent with their value of (47.7 ± 1.5)°. Es- 
sential for that agreement is, of course, not to consider 
the K — TTTT datum. Also, all the D wave threshold pa- 
rameters are now in good agreement with those used in 
[m . The remaining differences with respect to that work 
are rather small: the largest one is a 2.1 standard devia- 
tion disagreement, with respect to their predicted value 
6i = (5.67 ± 0.13) X 10-3. In general, and up to 500 
MeV, the results of [ll| fall within roughly one standard 
deviation of our analysis. For instance, at the kaon mass, 
our CFD SO wave phase shift is sI^\Mk) = 39.1 ± 0.6°, 
identical to theirs to the last digit, but our S2 wave is 
S^iMx) = -8.2 ± 0.6°, 0.3 degrees more than theirs, 
which is half an standard deviation. This good agree- 
ment does not deteriorate much above that energy. For 
instance, at 800 MeV, which is their matching point be- 
tween the calculated phase shifts and their input, they 



use an input value of 



(0) 



82.3 ± 3.4°. In contrast, 



we obtain 6^°^ = 85.2 ± 0.5° directly from the CFD set, 

whereas we find — 85.7 ± 1.6° when using the same 
CFD set inside GKPY equations, that is, one of their 
standard deviations. Above 800 MeV their amplitudes 
are part of the input and not solutions of Roy equations. 

Finally, we would like to remark that our best values 
for the scalar scattering lengths in Eq. ([^5]). are in very 
good agreement with the experimental results from pio- 
nic atoms [H, [13], that yield: 



^0 
.(0) 



,(2) 



0.280 ± 0.013(St.) ± 0.008(Syst.) M; 
— Of;4+0 033 A^-i 



or Xstt decays [38 



,(0) 



,(2) _ 



= 0.2571 ±0.0048(St.) 

±0.0025(Syst.) ± 0.0014(Ext.) 



Had we used them as additional constraints with the sta- 
tistical and systematic errors added linearly as we did 
with other decays, the difference with our best results 
would have been barely modified. 
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As we commented in Sect. IIII Bl the phase difference 
4°^(M|)-4^^(M|.) =_(52.5±0.8exp±2.8exp)° has been 
recently re- analyzed [2^. This is a considerable shift 
from the previous value of (57.27 ± 0.82oxp. ± 3rad. ± 
lchPTappr.)°j in much better agreement with ours and 
other previous dispersive analyses. Note that the new 
number is also in good agreement with our results in Ta- 
ble |IVl 



VIII. SUMMARY 

In this work, we have presented the derivation of a 
once-subtracted set of Roy-like dispersion relations — the 
GKPY equations. We have shown and explained that 
above 450 MeV, and up to 1115 MeV, they provide 
stronger constraints on tttt scattering amplitudes than 
other existing sets of dispersion relations. 

We have then applied these new equations as con- 
straints in our fits to data — together with the standard 
Roy equations and Forward Dispersion relations — in or- 
der to obtain a precise description of tttt scattering ampli- 
tudes. In contrast to previous works, we have extended 
the Roy and GKPY equations analysis from 932 MeV up 
to their applicability limit of 1100 MeV. Forward Disper- 
sion Relations are considered up to 1420 MeV. 

We have also made use of the final and very precise 
data on Ki4 decays from NA48/2, including the isospin 
violation corrections proposed in and we have re- 
moved a conflicting data point from K 2tt decay. 
With these changes in the data selection, most of the dis- 
agreement with previous Roy equation calculations [llj | 
has disappeared below 800 MeV. The largest discrepancy 
that remains is on the P wave slope parameter, but just 
at the two standard deviation level. 

In addition, we have improved our SO wave 
parametrization to ensure a continuous matching be- 
tween the low and intermediate energy parametrizations. 
Both parametrizations have been made more flexible, 
which allows the phase and inelasticity to include con- 
tributions from states different from tttt and KK, above 
the KK threshold. 

There are two sets of fits to data: unconstrained 
(UFD), or constrained with dispersion relations (CFD). 
In the UFD set each wave is independent of all others, but 
dispersion relations are satisfied only up to the two sigma 
level (in the sense explained in the text). In contrast, the 
CFD waves are all correlated, but they fulfill all disper- 
sion relations under consideration within less than one 
standard deviation in the whole energy region. The CFD 
set can be considered as a very precise parametrization 
of experimental data consistent with the requirements of 
analyticity, unitarity and crossing symmetry. Using this 
CFD set as an input in different sum rules and the disper- 
sion relations themselves, we have also provided a precise 
determination of phases in the elastic regime, threshold 
parameters and Adler zeros. 

In addition, and concerning the conflicting data for the 



SO wave inelasticity between the two-kaon threshold and 
1100 MeV, the use of the new GKPY equations has al- 
lowed us to show that the sudden drop around 1050 MeV 
in the SO wave inelasticity, or "dip solution", is clearly 
favored with respect to the "non-dip" solution. Actually, 
for the "non-dip" inelasticity scenario to fulfill dispersion 
relations, it would require a very poor description of the 
phase shift data, even when allowing for large systematic 
uncertainties. 

In conclusion, we provide fits to data in terms of sim- 
ple and ready to use parametrizations for the SO, S2, 
P, DO, D2 and F partial waves, between threshold and 
1420 MeV. Additional simple Regge parametrizations are 
given above that energy. In particular, the CFD set 
satisfies remarkably well all the analyticity and crossing 
symmetry constraints in the form of once and twice sub- 
tracted Roy equations and Forward Dispersion Relations. 
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Appendix A: Partial wave parametrizations 

In the following, we provide the parametrizations we 
use for each partial wave, and then, the parameters for 
the UFD and CFD sets. For brevity, we do not explain 
again why a specific parametrization for each wave has 
been chosen, since such details can be found in KPY08 
[1]. In what follows we use AU = 139.57 MeV, Mk = 496 
MeV and M„ = 547.51 MeV. 



1. SO wave 

This wave has been thoroughly discussed in the main 
text. However, for the sake of completeness, we repeat 
here the form of the parametrizations and provide the 
values of the parameters for the UFD and CFD sets in 
Table El 

For this wave we have set the matching point be- 
tween the intermediate and low energy parametrizations 
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clt — 0. 

use: 



cot4"^(s) 



85GeV. Thus, at low energies s < sm, we 



2fc 



2-^0 



X 



(Al) 



10(5) = 



Bo + Biw{s) + B2w{sf + B3w{sf 



So = 4M|.. 



(A2) 



Above that energy, and up to 1.42 GeV, we use the 
KPY06 polynomial parametrization for the phase shift, 
but with one more term in the expansion. For definite- 
ness, we provide here the polynomial parametrization 
once it has been matched to Eq. (jAip above, by impos- 
ing continuity and a continuous derivative at s = Sj\/, 
namely 



dn 1 



K 



+ |A:2|(fcA/- |fc2|) 



M 



K 



M2' 



^ {kM-MY ^ , (0.85 GeV)2 < s < AM^, 



Mi 



4M|. < s < (1.42GeV)2, 

m 



where ^2 = yJs/'i-Mf^. Note that we have defined 
5m = S{si\i) and (5^^ = d6{sM)/ds, which are obtained 
from Eq. (KT\ . and kM = |fc2(sA/)|. 

Finally, we assume an elastic SO wave, i]^^ = 1, up to 
the two-kaon threshold, whereas above that energy, we 
use: 



Vo (s) = exp 



-fc2(g) 

,sl/2 



-646 



4M, 



s 

2^,^3(5) 
J) j 



,1/2 



We have collected the values of the parameters for the 
UFD and CFD set in Table IVl 



SO wave 


UFD 


CFD 


Bo 


7.26 ±0.23 


7.14 ±0.23 


Bi 


-25.3 ±0.5 


-25.3 ±0.5 


B2 


-33.1 ± 1.2 


-33.2 ± 1.2 


Bs 


-26.6 ±2.3 


-26.2 ±2.3 


zo 






do 


(227.1 ± 1.3)° 


(226.5 ± 1.3)° 


c 


(-660 ± 290)° 


(-81 ±294)° 


B 


(94.0 ±2.3)° 


(93.3 ±2.3)° 


C 


(40.4 ±2.9)° 


(48.7 ± 2.9)° 


D 


(-86.9 ±4.0)° (-88.3 ±4.0)° 


?! 


4.7 ±0.2 


4.9 ±0.2 




-15.0 ±0.8 


-15.1 ±0.8 


h 


4.7 ± 2.6 


4.7 ± 2.6 


£4 


0.38 ±0.34 


0.32 ±0.34 



TABLE V: SO wave parameters for the UFD and CFD sets. 
The first four lines correspond to the low energy parametriza- 
tion, ^/s < 0.85 GeV , and the last nine to the parametrization 
up to ys = 1.42 GeV. 



2. S2 wave 

As we have already done with the SO wave, we have also 
set the matching point between intermediate and low en- 



1/2 



ergy parametrizations for this wave at s^^ 
Thus, at energies s^/^ < s^/j^ we use 



850 MeV . 



cot^^^^(s) 



wi{s) 



si/2 Mi 



2k s- 



2zl 



{Bo + B^wiis)} , 



,1/2 



/s + ^/si - s 



1.05 GeV, (A5) 



whereas at intermediate energies, 850 MeV < < 
1420 MeV , we use: 



cot4^^(s) = 



{Bho + Bhi[wh{s) - Wh{sM)] + Bh2[wh{s) - Wh{sM)f] , 
where 



Wh{s) 



y/s- \/sh - s 



,1/2 



Bho = -Bo ± BiWi{sm), Bhi = Bi 



1.42 GeV, 
dwi{s) 



dwh{s) 



Bhi — Bi- 



si \/sh - Sm ( y/sM + V^h- sm 



S = Sm 

2 



Sh \/si - Sm \ y/sM + y/si - SM 



(A6) 



Note that, with these definitions, both the parametriza- 
tion and its derivative are continuous at the matching 
point. 

Note that we have explicitly factorized the Adler zero 
at SA = 2z|. For the unconstrained fit, Z2 is fixed to the 
pion mass. As explained in the main text in Sect. IV Bl 
we then calculate the Adler zero position using Roy and 
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GKPY equations, and feed the weighted average into the 
constrained fit. This change is very smaU in terms of the 
total values and uncertainties of other quantities, but it is 
relevant in the differences when calculating the fulfillment 
of GKPY equations. 

For the inelasticity, we set 770 (s) = 1 for s below s — 
(1.05 GeV)^ and above that energy we use the empirical 
fit: 



P wave 


UFD 


CFD 


Bo 


1.055 ±0.011 


1.043 ±0.011 


Bi 


0.15 ±0.05 


0.19 ±0.05 


Ai 


1.57 ±0.18 


1.39 ±0.18 


A2 


-1.96 ±0.49 


-1.70 ±0.49 


ei 


0.10 ±0.06 


0.00 ± 0.06 


£2 


0.11 ±0.11 


0.07 ±0.11 



TABLE VII: P wave parameters for the UFD and CFD sets. 



The S2 wave parameters for UFD and CFD sets are 
given in Table IVII 



3/2 



S2 wave 


UFD 


CFD 


Bo 


-80.4 ± 2.8 


-79.4 ±2.8 


Bi 


-73.6 ± 10.5 


-63.0 ± 10.5 


Z2 




143.5 ± 3.2 MeV 


Bh2 


112 ±38 


32 ±38 


e 


0.28 ±0.12 


0.28 ±0.12 



TABLE VI: S2 wave parameters for the UFD and CFD sets. 



P wave 



For this wave we have set the matching point between 

1/2 

low and intermediate energy parametrizations at Sj(.j — 
2Mk- Thus, at low energies s^/^ < 2Mk, we use: 



cot 6i{s) 



„l/2 

2fc3^ P 



(2Ml_ 



Bo + Biw{s) \ , 



w{s) 



\fs - Vsq - s 



,1/2 



1.05 GeV, 



(A7) 



where the p mass is fixed to Mp = 773.6 ± 0.9 MeV. At 
intermediate energies, 2Mk < < 1420 MeV, we use 
a purely phenomenological parametrization: 

6iis) = Ao ± Ai {V^/2Mk - 1) ± A2 (^/i/2MK - l)' , 
mis) = l-ei^l-4M^/s-e2(l-4M^/s), (A8) 

where Aq is fixed from the value of 5i(4M^) obtained 
from the low energy parametrization, so that the phase 
shift is continuous. Note the possible presence of a dis- 
continuity in the derivative, allowed by the presence of 
the KK threshold. The values of the UFD and CFD 
parameters are given in Table IVlIl 



4. The DO wave 

As it was the case for the P wave, the matching en- 
ergy between low and intermediate energies is now taken 



1/2 



2Mk- At low energies, s^/^ < 2Mk, we 



at s 

parametrize this wave by: 
1/2 

cot 4°' is) = ^iK - s)Ml {Bo ± BMs)} : 



wis) 



/s ± - s' 



,1/2 



= 1.05 GeV, (A9) 



where the mass of the /2(1270) resonance is fixed at 
M/2 = 1275.4 MeV. In the intermediate region, 2Mk < 
gi/2 < 1420 MeV, we use a rather similar parametriza- 
tion: 



cot S. 



Whis) 



,1/2 



is) = ^iM^ - s)Mi {Boh + BihWhis)) ., 



y/sh 



,1/2 



1.45 GeV. (AlO) 



Imposing continuity at the matching point fixes Bho from 
the value of 4°^(4M|.) obtained fr om the low energy 
parametrization. We take the inelasticity to be differ- 
ent from 1 only for s > 4Af^, in which case we write: 



(0) 



l-e 



1 



1 -4M|./s 
4Mi/M2^ 



5/2 



k2is) 



(All) 

The parameters of the DO wave are given in Table IVIIII 



DO wave 


UFD 


CFD 


Bo 


12.47 ±0.12 


12.40 ±0.12 


Bi 


10.12 ±0.16 


10.06 ±0.16 


Bhi 


43.7 ± 1.8 


43.2 ± 1.8 


e 


0.284 ±0.030 


0.254 ± 0.030 


r 


2.54 ±0.31 


2.29 ±0.31 



TABLE VIII: DO wave parameters for the UFD and CFD sets. 



5. The D2 wave 

We use the following parametrization from threshold 
up to 1420 MeV : 

cot4'^(s) 



2A;5 



wis) 



4(M2±A2)-s 
_ - y/sp - s 

y/s + - S ' 



{Bo + Biwis) + B2wisf} , 



,1/2 



1.45 GeV, (A12) 
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and we consider that the inelasticity differs from 1 for 
> 1.05 GeV, as follows: 

ryf'(s) = l-e(l- J/s)3, s^/^ ^ 1.05 GeV , (A13) 

which is almost negligible up to 1.25 GeV. The values 
of the parameters for the UFD and CFD sets are given 
in Table HXl 



D2 wave 


UFD 


CFD 


Bo 


(2.4 ±0.5) 10^ 


(4.1 ±0.5) 10^ 


Bi 


(7.8 ± 1.0) 10^ 


(8.6 ± 1.0) 10^ 


B2 


(23.7 ±4.2) 10^ 


(25.5 ±4.2) 10^ 


A 


196 ± 25 MeV 


233 ± 25 MeV 


e 


0.2 ±0.2 


0.0 ±0.2 



TABLE IX: D2 wave parameters for the UFD and CFD sets. 



6. The F wave 

We neglect the inelasticity up to 1420 MeV and simply 
use the following parametrization from threshold: 

cotS,is) = +Bo + BM^s)^ , 

v^-7i^ = 1.45 GeV. (A14) 

The parameters for the UFD and CFD sets are given in 
Table [Xj Note that they do not change at all from one 
set to another. 

F wave UFD CFD 

Wo (1.09 ± 0.03) lO'' (1.09 ± 0.03) lO'' 
Bi (1.41 ±0.04) 10^ (1.41 ±0.04) 10^ 
A 0.051 X 10^ 0.051 X lO'^ 



TABLE X: F wave parameters for the UFD and CFD sets. 



Regge param. 


UFD 


CFD 




1.22 ±0.14 


1.48 ±0.14 


MO) 


0.46 ± 0.02 


0.53 ±0.02 


Pp 


2.54 ± 0.04 


2.50 ± 0.04 


Cp 


0.0 ± 1.0 GeV -2 


0.6 ± 1.0 GeV -2 


Cp, 


-0.4 ±0.4 GeV 


-0.38 ±0.4 GeV 


pp' 


0.83 ± 0.05 


0.80 ± 0.05 


ap,{Q) 


0.54 ± 0.02 


0.53 ±0.02 




0.2 ±0.2 


0.08 ±0.2 


Bp 


0± 2.5 GeV-'' 


2.7 ±2.5 GeV-* 



TABLE XI: UFD and CFD parameters for the p, Pomeron 
and 1 = 2 Regge contributions to tttv scattering amplitudes. 



4Af^ > t > — 0.4 GeV ^, somewhat less rehably for the 
most negative t values. This is enough to describe the 
region of interest, that reaches t — —0.42 GeV^. In par- 
ticular, for the p Regge trajectory, we use the following 
expression for the imaginary part, which is all we need 
in the dispersive integrals: 

ImF(^*-)(M) = /3.^±^^We-(f)"^^r 
l + ap(0) \s/ 

ap{t) = ap{Q) ±iap ± ^t^a", 

if{t) = l + dpt + ept^, (A15) 

where we fix: 

s = lGeV^ 6 = 2.4 ±0.2 GeV (A16) 
a'p = 0.90 GeV = -0.3 GeV -^ 
dp = 2.4 ±0.5 GeV -^ 

whereas the rest of the parameters are allowed to vary in 
the fits. 

For both the Pomeron P and the P' pole, we have used 
for = 1420 MeV : 

ImF(^*=°)(s,i) = P{s,t) + P\s,t), 



7. The G waves 

The contribution of the GO and G2 waves was shown to 
be completely negligible for the calculations. The details 
can be found in the Appendix of KPY08 [3|] . 



8. Regge Parametrizations 

Next we show the Regge parametrizations that we use 
in the high energy region, i.e. above 1420 MeV . The for- 
ward {t = 0) Regge parametrizations were obtained from 
fits to high energy data [1]. For the t ^ behavior we 
simply covered the uncertainties between the different 
fits in [3J]. These Regge fits are expected to represent 
experimental data when 1.42 GeV < s^/^ < 20 GeV and 



Pis,t)^l3p'<SJp{t)apit) 
ap{t) — 1 + ta'p, 



l + ap(t) 



1 



Xp{t) 



^p{t) = 1 ± Cpt, 



ap,{t)[l + ap,{t)] 



ap,{t) 



P'{s,t) = l3p,^p,{t) --— --.e' 

Q;p'(0)[l + ap'({))\ \s 

ap,{t) = ap,{0) +ta'p,, fp^i) = I + cp,t, 

(A17) 

where, once again, we fix: 

s = lGeV^ 5 = 2.4 ±0.2 GeV (AlB) 
a'p = 0.20 ±0.10 GeV -^ a'p, = 0.90GeV-^ 
Cp = 0.0 ± 1.0 GeV -^cp, = -0.4 ±0.4 GeV 

and allow the rest of the parameters to vary in the fits. 



29 



Finally, the Regge exchange of isospin two is 
parametrized as: 



ImF 



s\"p(«)+"p(o)-i 



(A19) 



In Table I XII we show the values of the Regge param- 
eters obtained from the direct fit to high energy data 
(UFD) and how they are modified when imposing the 
dispersive constraints in the fits (CFD). 



Appendix B: Derivation of the once subtracted 
dispersion relations 

A once subtracted dispersion relation for a scattering 
amplitude of definite isospin / has the following expres- 
sion: 



S - So 



s-so Im F('Hs',t) 

4M2 {s' - sa){s' - s) 

, Im i^(^)(s',i) 



ds' 



(s'-so)(s'-s) 



, (Bl) 



with So the subtraction point. This expression assumes 
that the point s is regular. However, we are especially 
interested in what happens for s in the physical region, 
that is, on the cuts of the function F{s,t). The usual 
prescription is to define the amplitude for physical values 
of s as: 



Fphys{s,t) = lim F{s + ie,t). 

e->0+ 



With this prescription, we have: 



s — So + ie 



lmF'-^Us',t) 
as — - 

4M2 (s' - so)(s' - s - ie) 

s - so + ie r°°^^, ImF(^)(s',i) 

TT J-t 



(s' - So)(s' - s~ie)' 



To obtain the physical amplitude, we must take the limit 
e — > 0+ in this expression. Suppose s is on the right hand 
cut (RHC), 4M2 < s < oo. Since 



1 



= P.P. 



=F iTrS{x), 



0+, 



we can write the RHC integral as: 

s-so f°° , lm.F^^\s',t) ,r., , 
^ P.R / ds' — ..\ \ + I Im F(^) s, , 



whereas the left hand cut (LHC) integral presents no 
problems when e vanishes. Then we have 

F'£u,t)^F^'\s,,t)+^l^ F('Hs,t) 



+ 



So 



P.P. 



oo 

S - So 
TT 



ds 



, IniFW(s',t) 
(s'-so)(s'-s) 

, Im i^(-f)(s',t) 



ds' 



is' - so)is' - s)-' 



thus the dispersive integrals only reconstruct the real part 
of the amplitude, instead of the total amplitude. Had 
we chosen s to be on the LHC, the reasoning would be 
analogous, but the principal value should be taken on 
the LHC integral, instead of on the RHC one. We finally 
obtain: 



ReFl,il{s,t)^ReF('Hso,t) 



s - So 



P.P. 



ds 



s - So 



4M2 



P.P. 



, Im F^''>{s',t) 
(s'-so)(s'-s) 

, Imi^W(s',t) 



ds' 



is' ~ So)is' - s)-' 



with the principal value taken on the cut on which s lies. 
This is valid for any s on the cuts of F^-^^s^t), i. e., for 
physical s. We can now recast the LHC integral on the 
s-channel in Eq. (IB1[) in terms of the u-channel RHC by 
renaming the dummy variable s' as u' in the LHC integral 
and performing the substitution 

u' AMI -s' -t. 

Taking both integrands under the same integral sign, and 
choosing so — — in analogy with Roy's derivation — we 
obtain: 

Re F^'\s,t)=Yie F^'\o,t) 



ds' 



TI" J AMI 



linF^^\s',t) lraF^^\u',t) 
s'{s' — s) u'{u' — s) 



Each of these integrals is potentially divergent if taken by 
themselves due to the Pomeron contribution coming from 
the It = channel, which grows like Im F(^'=")(s,t) ~ s 
for large s. We now show that this is not the case when 
taken together. 

Bose statistics require that the It — amplitude be 
symmetric under s — u exchange, 

Since the amplitudes with well-defined isospin in the s- 
and i-channels are related via the usual crossing matrices, 

/l/3 1 5/3 \ / 1/3 -1 5/3^ 

Cst = 1/3 1/2 -5/6 , Csu = -1/3 1/2 5/6 
\l/3 -1/2 1/6 J \ 1/3 1/2 l/6y 
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we know that each amphtude with well-defined isospin 
in the s-channel has a contribution from each of the am- 
plitudes with well-defined isospin in the i-channel. In 
particular, the contribution from the It ^ channel to 
the integrand can be written as: 



1 



1 



s'{s' — s) u'{u' — s) 

{s + t- 4M2)(2s' + t - 4M2)Im F'^^^^^^s' ,t) 
s'is' - s){s' + t ~ AM^){s' + s + t- AMI) ' 



The s'^ terms in the numerator cancel out, and the in- 
tegrand decays as when s' ^ oo, so that the in- 
tegral converges. This is in contrast with the expected 
1/s' asymptotic behavior, which would spoil convergence. 
The contributions from the other i-channel isospin con- 
tributions It = 1,2 are not problematic, since they grow 
as (s')" with a < 1, and are convergent even if taking the 
integrals separately. Note that this cancellation does not 
depend on the explicit parametrizations we use for the 
Pomeron, but rather, on very general asymptotic prop- 
erties of the amplitudes. 

In order to rewrite the RHC contribution from the u 



channel in terms of amplitudes on the RHC s-channel, 
we take into account the crossing symmetry relation: 

f(^)(4M2 - s' - t,t) = Y,Cii'F^''\s',t), (B2) 

with Csu the crossing matrix defined above. Also, 

F<-'H0,t)^J2^""F'^'"Ht,0): (B3) 
/" 

and we now write a dispersion relation for F^-^ \t, 0): 



F(^")(i,0)=F(^")(to,0) 



t-to 



ds' 



Im i^(^")(s',0) 
[s' - t){s' - t^) 



[AMI ~ t - s'){AMi ~ s' - to) 
Again, in analogy with Roy, we take to = AM^. Thus: 



(B4) 



poo 

Re F^'\s,t) = Y Cl[ F(''\AMI,Q) + -RR / ds' 



lmF('Hs',t) ErCf^ImF(^')(s^^) 
s'{s' -s) is' + t^ 4M2)(s' + s + t- AM^) 

Im F(^")(s',0) E/'" Cf:^"'lm F('"'\s',0) 



[s' - t){s' - AMI) 



s'(s' + t-4M2) 



Now, to project into partial waves, we define first the where Pi{x) and Pe'iy) are Legendre polynomials, and 
following kernels: 

(,s-4M2)(x-l) 

s * = 9 ' 

Ku'{s,s')=— dxPe{x)Pe> (y) , (B5) , 

u' = AM'^-s'-t, 



y 



Lu'{s,s') = ^ I dxP,{x) ^''^y} " 

^ '^l Note we have taken advantage of the symmetry of the 







1 t — AM^ integrands to change the integration limits from (—1, 1) 

^^(^'^)^ 4^--4M2) X ^"^^(")^73F' to (0,1). 

_^ I t With the normalization chosen in Section III Bl and re- 

Ni{s,s') = / dxPi{x) , calling that al^^ = 0, we find: 

TTS' Jo u' 



poo ( ^ 

Re t^p{s) ^ ^iV] CU"a[' ^ + V(2£' + 1) / ds' \ Ku'{s, s')Im tf(s') ~ Lu'{s, s') V C.^Jm t[! \s') 
7^ V •^'iA/^ [ V 



M,(s, s')Im tf\s') ~ mis, s') C^I"/'"Ini 4^"'\s') 
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In order to simplify the previous expression, we define 

(B6) 



Ku' (s, s') 



We thus arrive at the final result used in Eq. ([18 



= (2/ + 1) 
X [Ku,{s,s')S''' -Lu'{s,s'){Csuy'' 
+ Me{s, s'){Cl[) - iV,(s, s'){CstCsuY''] 



Re t[^\s) 



= ST,+DT,{s) 



4M2 



ds'i^„, (s,s')Im 4'{s'), 



where, for simplicity, the high energy part of the inte- 



grals (s' > Smax) and the higher partial waves {£' > 1) 
are grouped in the so-called driving terms, DT^[s). The 
subtraction terms ST^ , which are now constant, are: 



with the coefficients defined in Eq. (|B7p . For our pur- 
poses we will only need Co = 1 and — 1/2. Note that 

the subtraction term STg is a constant, and does not de- 
pend on s. This is a relevant feature of GKPY equations 
versus Roy Equations, as explained in Section [IV El 



dxPAx) = — - — 

2r(i 



mm 



1, 

0, 



(-1)" 



£^0 

£ — 2to, m > 



2™+i(m+l) 



TV UT=o [2™ - (2fc + 1)] , £ = 2m+l 



(B7) 



Appendix C: Integral kernels in GKPY equations 



All kernels in Eqs. (IB5P - (|B6p can be calculated ana- 
lytically. One has to note, however, that the L«'(s,s') 
and Ni^Sjs') kernels are singular at u' = 0, namely 
X = — (2s' — s — 4M^)/(s — 4M^), where a principal value 
over the integral is understood. 

In this work we need 18 Kg^, (s, s') kernels, since we 
are considering the dispersion relation for the SO, P and 
S2 waves, but using SO, P, S2, DO, D2 and F waves as 
input. However, following (soj . we know that, since the 
K, L, M and N kernels in Eqs. ([B5])-([B6| do not depend 

on isospin, the K^/f,, (s, s') are not all independent and 
can be expressed in terms of four of the Kw above, and 
eight combinations of the other kernels, which we call 
/«'(s,s')- Namely: 



-"■00 



-"■01 



7f2 



3^01 , 



/oo/3. 



nn 



-3/00, 



25 
-y/02, 



kII = 7/03, 



i^io 



/10/3, 

3(Xii 

25 
-y/12, 



K 
K 



-g/10, 

I,. 



kII = 7(Xi3 - i/13), 



Kqo — -loo/ 3, 
-21 3 



-^00 ^ -^00 — ^00/6, 



Iqi, 



2 

5(i^02 



■20 
02 

■21 

03 



02, 



where 



loo = Lqq 
ho = LiQ 
I02 — L02 



Mo- 
Ml 



No, 


hi 


— ioi 


+ Mo 


-No, 


Ni, 


hi 


= Lu 


- Ml 


-Ni, 


No, 


I03 


= -^03 


+ Mo 


-No, 


Ni, 


/l3 


= L13 


- Ml 


^Ni. 



The analytic expressions for the Kn^ kernels are: 

Koo ' 
K02 

Ki3 



ns'{s — s') ' 
s(4M2 + .s - 2s') 
" 27rs'(s'-4M2)2 



s(8M2 



3s') 



67rs'(s-s')(s'-4M2)' 



s(4M2 



2s' 



r\2 



87rs'(s' - 4M2)3 ' 



(CI) 



(C2) 



The diagonal kernels Koois,s') and Kii{s,s') contain a 
singularity at s = s', which is the only type of singularity 
in the GKPY equations. 
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By defining the following functions: the analytical expressions for the Iu'{s, s') can be recast 

as: 

= s + s'-AMr ^^'^ 
(s + 2s' - 4M2)2 

= Ms + s'-4M^)^ ' J _ o is - ^Ml){s' - 2Ml)/{s' - AMI) + s' log (ai) 

^ - 4{s' - 2Mlf irs'is-AMl) 

4(s'-4M2)(. + ,s'-4M2)' ^^^(,^,/) = .^(flzl^ 

(s - 2s' + 4Af2)(,s + ,s' - AMI) ~ 



«4 - (s'-4M2)(s + 2s'-4M2) ' o {s' - ^Ml)s' \og (ai) + ss' \og (a2) ^ 

as 



Ct6 



+ 2s' - AMI) - - 



(s'-4M2)(s + 2s'-4M2)' 
(s - 2s' + AMl)(s + 25^ - 4M2) 
4(s' - 4M2)s' 



r / A 1 / 6s 1 1\ 1 /„, , ^ 6s(s + s'-4M2)log(a2)\ 

1 f (s - 4M2) (2s'3 + 10(s - 2M2)s'2 + (25s2 - SOM^ s + 64M^) s' - 64M^) 



-fo3(s,s') 



7r(s - 4M2) 1 (s' - 4M2)3s' 

^ , 2s (I0s2 + 15 {s' - 4M2)s + &{s' - AM^)^) log ja^) \ 
21og(ai) + (a' - 4M2)3 1' (C^) 



-fio(s,s') 



7r(s-4M2)2(s'-4M2)s' 

+ (s' - 4M2)s' log (ai) s + 2s' (s'^ - Qs' Ml + 8M^) log (ai) 



s^Ml + 2s'^s - Ss'sM^ - 8sM^ - Ss'^M^ + 32s'M^ + 16M^ 



(C8) 



(2s' + M2)s2 + 2 (s'2 - 8 s'M2 - AM*) s - 8 (s'^ - 4M2s' - 2M^) M^ 



7r(s - 4M2)2 
1 



(s' - 4M2)s' 

s(s + 3s' - 8Ml) log (as) - (s^ + 2 (s' - 2M2)s + 2 (s'^ - 6s'M2 + 8M^)) log (ai) 



s'-4M2 

+ 2(s'-2M2)log(a5)-slog(a6) 



(C9) 



/12(S,S') = 



1 



2(s - 4M2) (9s's2 + 2 (6s'2 - 17 s'M2 - 4M^) s + 4 (s'^ - 8s'2m2 + 14M^ s' + 8M^)) 



27r(s - 4M2)2 
4 

(s' - 4M2)2 . 

(3s3 + 3(3s' - 8M2)s2 + 6 (s'^ - 6s'M2 + 8M^) s + 2(s' - 4M2)2 (s' - 2M2)) log (04) 



(s' - 4M2)2s' 

s (3s2 + 3(3s' - 8M^)s + 7s'^ - 44M^s' + 64M^) log (a^^) 



- 8(s'-2M2)log(a5)+4slog(a6) }, 



(CIO) 
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' 2(s - AMI) (85s's3 + 5s'(33s' - lOOAf^) + (J2s'^ - SlOs'^M^ + IMM^s' + 96Af6) s) 

3(s' - 4A//2)3s' 

lOs^ + 5(7s' - 20AfJ) + 12 (Ss'^ - lOs'Af^ + 28M^) s 



(s' - 4Af2)3 

{s' -AMl)\Us' ~2SMl) \ log( 



a. 



+ ( 10/ + 5(7 s' - 2{)Ml)s^ + 12 (3s'2 - 19s' A/^ + 28Af^) + 12(s' - AMI)'^{s - 2M^)s 



+ 2{s'~AM^)^{s'-2M^)j\og{a,) 
s - AMI) (s" - 4A^^s' _ 2Af4) 



+ 



(s' - 4Af 2)s' 



+ 8(s'-2Af2)log(a5)-4slog(a6) 



(Cll) 



The behavior around threshold is also interesting when 
considering the expansions of the kernels around s — AM^. 

In particular, the threshold expansions of Kgg, (s, s') 
around s = 4Af2 behave like a + b{s - AMI) + • ■ • 



do not weight the uncertanty but take the smallest of the 
two outputs, since both results come from the same data. 



Appendix D: ROY-GKPY weighted phases 

In Table IXIII we give the central values of the phase 
in the elastic regions, as the weighted averaged obtained 
from the output oj Roy and GKPY equations, when using 
the CFD set as input. These results could be understood 
as a traditional "energy dependent data analysis". We 
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940 
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970 
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-23.9±1.7 



TABLE XII: Phases from the dispersive data analysis. Cen- 
tral values are obtained as a weighted average between the 
output of Roy and GKPY equations, using the CFD fit as in- 
put. We do not weight the uncertainty but take the smallest 
of the two, since both results come from the same data. 
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